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FOREWORD 


This  document  is  the  Master's  Thesis  of  Roger  C.  Hendershot.  It  is 
being  issued  as  a  technical  report  for  the  Advanced  Lightweight  Torpedo 
project  under  Contract  N00024-78-C-6018.  The  work  reported  is  part  of 
a  study  of  digital  simulations  of  high-frequency  acoustic  backscattering 
characteristics  of  fish  schools  and  bottom  features.  The  aim  of  this 
study  is  to  provide  accurate  simulation  of  false  target  echoes  in  the 
Applied  Physics  Laboratory's  Reverberation  Generator  (REVGEN)  simulator 
and  the  Naval  Ocean  System  Center's  Hybrid  Simulator. 

This  study  is  in  an  early  stage,  and  the  work  reported  here  relates 
to  a  relatively  simple  fish  school  model  designed  to  aid  in  the  inter¬ 
pretation  of  field  data  and  to  provide  insight  into  some  of  the  problems 
of  digital  simulation  of  fish  school  echoes.  As  the  project  continues, 
detailed  point-scatterer  models  for  biological  and  bottom-feature  false 
targets  are  to  be  developed,  realized  in  software,  and  tested  by  com¬ 
parison  with  field  data. 

Progress  made  on  a  detailed  fish  school  model  will  be  documented  in 
another  APL-UW  report,  "Construction  of  a  Simulated  Fish  School  as  a  Sonar 
False  Target,"  by  G.E.  Lord.  That  document  contains  a  description  of  the 
model  and  a  listing  of  the  REVGEN  software  that  incorporates  the  model. 

In  that  model,  each  scatterer  follows  a  separate  track  through  the  ocean, 
each  track  having  both  deterministic  and  random  components.  The  Doppler 
shift  and  aspect-dependent  cross  section  are  computed  from  tracking  infor¬ 
mation  . 

In  contrast  to  the  REVGEN  simulation,  the  model  described  in  this 
thesis  is  simpler  in  that  there  is  no  tracking  or  "memory"  of  fish  po¬ 
sitions  on  successive  pings.  This  simpler  model  serves  two  purposes: 

(1)  it  is  a  convenient  means  for  determining  the  sensitivity  of  the 
simulation  output  to  parameters  such  as  number  of  point  scatterers, 
distribution  of  individual  scattering  cross  sections,  distribution  of 
Doppler  shifts,  etc.,  and  (2)  it  provides  a  method  of  determining  biases 
in  fish  school  size  estimates  made  using  acoustic  data.  An  understanding 
of  these  biases  is  important  in  extracting  depth  and  size  distributions 
for  the  fish  schools  in  the  experimental  data  base. 

This  thesis  emphasizes  split-beam  processing  methods,  more  partic¬ 
ularly,  crosscorrelations  between  split  beams.  This  choice  was  made 
because  a  rather  large  quantity  of  field  data  has  already  been  processed 
using  this  technique,  with  promising  results.  Also,  we  have  found  that 
the  crosscorrelator  output  is  extremely  sensitive  to  changes  in  simula¬ 
tion  parameters  and  provides  a  clear  indication  of  the  onset  of  trouble 
when  scatterer  density,  sampling  rate,  etc.,  take  on  unacceptable  values. 


ABSTRACT 


A  computer  simulation  based  on  a  point  scatterer  model  has  been 
used  to  investigate  the  problem  of  estimating  the  location  and  dimen¬ 
sions  of  fish  schools  using  split-beam  processing  of  hydroacoustic 
data.  Estimates  are  derived  from  the  temporal  crosscorrelation  of 
signals  received  by  transducer  halves,  the  peak  correlation  amplitude 
being  related  to  the  variance  of  the  scatterer  locations  (the  fish 
school's  width)  and  the  time  shift  of  the  peak  being  related  to  the 
mean  of  the  scatterer  locations  (the  fish  school's  angular  position). 

These  estimates  are  formed  with  vertically  and  horizontally  separated 
transducer  halves  and  combined  with  echo  arrival  time  to  provide  three- 
dimensional  position  and  size. 

It  is  demonstrated  that  fluctuations  in  the  estimated  angular  size 
decrease  as  the  ratio  of  pulse  length  to  echo  length  decreases. 

The  effect  of  additive  noise,  uncorrelated  between  acoustic  channe.  5, 
has  also  been  studied.  Noise  causes  width  estimates  to  be  biased  toward 
larger  values,  owing  to  decorrelation  between  the  acoustic  channels.  The 
location  estimates,  however,  remain  unbiased  for  signal -to-noise  ratios 
down  to  approximately  -5  dB. 
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CHAPTER  1 


INTRODUCTION 


In  hydroacoustic  assessment  of  fish  populations,  the  ability 

to  retrieve  spatial  characteristics  from  fish  school  echoes  is 

highly  desirable,  buc  generally  difficult.  Estimations  of  the 

overall  school  size  and  location  are  of  interest  in  determining, 

for  example,  behavior  patterns  and  mean  target  strengths  of  fish  in 

their  natural  environment.  Narrowbeam  sector- scanning  arrays  can 

be  used  to  determine  spatial  properties.  The  use  of  split-beam 

systems  is  another  attractive  possibility  because  of  their  relative 

low  cost  and  simplicity  compared  with  sector-scanning  systems. 

Split-beam  systems  have  been  used  to  determine  school  depth1  and 

2 

target  strengths  of  single  fish. 

In  this  thesis,  computer  simulations  are  used  to  investigate 
split-beam  techniques  for  estimating  overall  school  size  and  loca¬ 
tion.  Chapter  2  provides  a  general  description  of  the  split-beam 
method.  Fish  schools  are  simulated  using  ellipsoids  filled  with 
point  scatterers,  and  echoes  from  these  schools  are  constructed 
using  the  single  scattering  approximation.  The  crosscorrelation 
of  these  echoes  received  by  the  quadrants  of  the  split-beam 
transducer  is  then  used  to  calculate  the  size  and  location 
estimates . 

A  description  of  the  simulation  programs  is  given  in  Chap¬ 
ter  3.  A  series  of  programs  is  used  to  study  the  effect  of  the 
school  characteristics,  noise,  and  sampling  parameters  on  the 
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spatial  estimates.  In  Chapter  4  the  results  of  the  simulations  are 
given  along  with  supporting  analysis  for  simplified  situations. 
Biases  in  the  dimension  estimates  resulting  from  the  decorrelating 
effect  of  the  noise,  and  other  sources  of  fluctuations  in  the 
estimates  are  illustrated.  Conclusions  about  the  errors  in  split- 
beam  spatial  estimates  follow  in  Chapter  5. 


CHAPTER  2 


DESCRIPTION  OF  THE  SPLIT-BEAM  SIMULATION 


An  introduction  to  the  split-beam  method  and  a  general  descrip¬ 
tion  of  the  computer  simulation  are  presented  in  this  chapter.  The 
formulation  of  crosscorrelations  of  the  echo  signals  using  quadra¬ 
ture  sampling  is  given  along  with  the  relationship  between  these 
crosscorrelations  and  the  angular  location  and  dimension  estimates. 

A  physical  representation  of  the  split-beam  technique  used  in 
the  computer  simulation  is  shown  in  Fig.  1.  This  system  uses  a 
multiple-element  transducer  with  the  elements  combined  to  form  the 
four  transducer-half  signals  labeled  left,  right,  up,  and  down.  In 
general,  the  target  will  not  lie  on  the  axis  of  the  transducer  beam 
and  the  four  signals  will  be  displaced  in  time  from  one  another. 

This  time  difference,  measured  by  using  the  time  shift  at  the  peak 
correlation  between  the  signals,  can  be  used  to  estimate  the 
bearing  to  the  center  of  the  target.  The  crosscorrelation  of  the 
left-right  signal  gives  the  horizontal  bearing  and  that  of  the  up- 
down  signal  gives  the  vertical  bearing. 

The  magnitudes  of  the  crosscorrelations  can  be  used  to  esti¬ 
mate  the  target's  angular  size.  The  greater  the  angle  covered  by 
the  target,  the  lower  the  crosscorrelation.  This  concept  can  be 
illustrated  with  a  simple  example.  A  single  scatterer  will  result 
in  identical,  time  displaced  signals  whose  normalized  -peak  cross¬ 
correlation  will  equal  unity-.  When  a  second  scatterer  is  added, 
the  arrival  time  of  its  echo  relative  to  the  other  scatterer's  will 


4 


differ  for  each  signal  as  can  be  seen  by  examining  the  path  lengths 
in  Fig.  2.  The  signal  envelopes  at  transducers  A  and  B  are  s^Ct) 
and  Sg(t),  where  Sg(t)  illustrates  the  possibility  of  destructive 
interference  between  the  two  scatterer  returns.  Since  these  sig¬ 
nals  differ,  their  normalized  peak  crosscorrelation  will  be  less 
than  unity.  Additional  scatterers  will  result  in  a  further  de¬ 
crease  in  the  peak  crosscorrelation. 

The  fish  school  is  represented  in  the  simulation  by  an  ellip- 

\ 

soid  filled  with  point  scatterers  whose  locations  are  randomly 
chosen.  In  the  simulations  reported  here,  the  ellipse  rotation 
angle,  0,  equals  zero,  so  the  x'  axis  will  always  point  to  the  cen¬ 
ter  of  the  transducer.  The  angle  ip  measured  in  the  left-right 
plane  of  the  transducer  is  the  target  bearing  relative  to  the  beam 
axis. 

The  number  of  point  scatterers  per  cubic  meter  in  the  ellip¬ 
soid  is  typically  around  0.3.  Multiple  scattering  effects  should 

3 

be  minimal  for  actual  fish  schools  with  these  densities.  The 
onset  of  first-order  multiple  scattering  occurs  when  the  term 
p<at>L  ~  0.1.  This  term  is  known  as  the  optical  distance  within 
the  scattering  medium,  where  p  is  the  density  of  the  scatterers, 
<at>  is  the  average  extinction  cross  section  of  the  individual 
scatterers,  and  L  is  the  rounl  trip  acoustical  path  length  in  the 
scattering  medium.  A  typical  optical  distance  simulated  in  this 
study  is  for  an  ellipsoid  with  x’,  y',  and  z  dimensions  of  10  m. 
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Figure  1.  Physical  representation  of  the  simulation. 
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Figure  2.  Echo  crosscorrelation  using  two  scatterers. 
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15  m,  and  5  m,  and  containing  100  fish.  For  this  situation  p  - 
100/ (4/3*Tr,2 ,5*7 .5*5)  =  0.255  fish/m3,  and  L  =  20  m.  The  extinc¬ 
tion  cross  section  of  a  fish  of  length  30- cm  is  approximately 

2  2 
30  cm  with  fluctuations  from  3  to  300  cm  depending  on  the  fish's 

orientation.  The  resulting  optical  distance  is  0.015,  which  is 

well  below  the  onset  of  first-order  multiple  scattering. 

The  point  scatterers  are  assigned  scattering  amplitudes  and 
Doppler  shifts.  The  scattering  amplitudes  either  are  equal  or  are 

4 

randomly  chosen  from  a  Rayleigh  distribution  using  a  different 
distribution  mean  for  each  scatterer.  The  resulting  amplitudes 
fluctuate  approximately  6  dB.  The  scattering  amplitude  of  each 
scatterer  is  also  corrected  for  the  strength  of  the  transmit  and 
receive  beam  patterns.  The  Doppler  shifts  are  either  zero  or 
chosen  from  a  uniform  distribution  ranging  from  20  to  100  Hz  using 
a  30  kHz  carrier  frequency  reflecting  relative  motions  ranging  from 
1  to  5  kn.3 

Four  echo  delays  for  each  scatterer  are  computed  using  path 
lengths  measured  from  the  center  of  the  transducer  to  the  scatterer 
and  back  to  the  left,  right,  up,  and  down  centers.  No  multiple 
scattering  paths  are  included.  The  effects  of  differential  spher¬ 
ical  spreading  and  absorption  losses  associated  with  the  path  lengths 
are  negligible  for  the  geometries  simulated  and  are  not  included. 

Next,  samples  of  the  fish  school  echo  are  constructed  as  shown 
in  Fig.  3.  Because  the  echo  signal  is  assumed  to  be  narrowband. 
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and  to  retain  phase  information,  samples  of  the  signals  are  taken 
in  quadrature  form.  If  the  echo  from  a  scatter er  is  present  at  the 
time  a  sample  is  taken,  the  scatterer's  contributions  to  the  quad¬ 
rature  components  of  that  sample  are  included  in  the  following 
summations. 

The  ith  quadrature  samples  taken  at  time  t^  are  given  by 
zL(i)  " 

E  A(k)BLR(k)P[t.-TL(k)]  exp{j2Tr[f0+fd(k)]-TL(k)( 

k 

zR(i)  - 

S  A(P,5BLR(J'.)P[ti-TRC?-j]  exp(j2lt[f0.fd(«j]-TR(?,)(  (1) 

zU(i)  = 

Z  A (m)  ByD (m) P  [t^-Ty (m) ]  exp  { j27r[f0+fd(m)  ]  ♦TU(m)  [ 
zD(i)  = 

E  A(n)BUD(n)P[t.-TD(n)]  exp  {  j27r[f0+fd(n) ] -cD(n) (  , 

where  the  summation  is  over  the  scatterer  echoes  k,  i,  m,  and  n, 
present  at  the  left,  right,  up,  and  down  transducer  centers  at  time 
t^;  A(  )  is  the  scattering  amplitude;  B(  )  is  the  combined  transmit 
and  receive  beam  correction  for  the  left-right  or  up-down  trans¬ 
ducer  halves;  P(  )  is  the  amplitude  of  the  scatterer  echo  envelope 
at  time  t^;  fQ  is  the  carrier  frequency  (30  kHz) ;  f ,  is  the  Doppler 
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shift;  t(  )  is  the  round-trip  delay  time  to  the  left,  right,  up,  or 
down  transducer  center. 

The  envelopes  of  the  scatterer  echoes,  denoted  by  P,  are 
assumed  to  be  identical  to  the  envelope  of  the  transmit  pulse,  and 
no  changes  caused  by  the  scatterer  or  the  medium  are  included. 

Also  note  that  because  of  the  difference  in  arrival  times  at  the 
four  quadrants,  a  scatterer  that  is  included  in  the  summation  for 
the  left  quadrature  sample,  for  example,  may  not  be  included  in  the 
summation  for  the  right  quadrature  sample;  hence,  the  different 
subscripts  k,  &,  m,  and  n. 

Uncorrelated  noise  is  included  by  adding  two  random  white 
Gaussian  noise  components  to  the  echo  samples.  The  signal-to-noise 
ratio  (SNR)  is  defined  for  each  echo  as  the  ratio  of  the  energy  of 
the  signal  samples  to  the  energy  of  the  noise  samples.  These 
signal-plus-noise  samples  are  then  digitally  filtered  by  a  single¬ 
pole  filter  whose  3  dB  bandwidth  is  nominally  1  kHz. 

The  four  sets  of  filtered  samples  are  each  normalized  by  the 
square  root  of  their  respective  total  echo  energies.  The  quadra¬ 
ture  samples  in  this  form  are  shown  as  Z L,  Z^,  Z^,  and  Zp  in 
Fig.  3. 

The  next  step  in  the  signal  processing  procedure  is  the  echo 
crosscorrelation.  In  general  the  crosscorrelation  of  two  sets  of 
echoes  would  be  a  function  of  the  shift  in  samples  between  the  two 
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sets.  However,  the  information  required  to  estimate  the  target's 
angular  position  and  dimensions  is  contained  in  an  expression  where 
no  shift  is  used.  This  expression  is 


CLr(°)  =  2  ZL(i)ZR*(i) 

cUD(°)  =  l  zu(i)zD*(i) , 


(2) 


x 

where  the  summation  is  over  all  echo  samples.  The  need  for  com¬ 
puting  only  one  crosscorrelation  term,  which  is  an  advantage  of  the 
split-beam  system,  can  be  understood  by  considering  the  crosscorre¬ 
lation  of  two  signals  c(t)  computed  in  continuous  time.  This 
crosscorrelation  is  related  to  the  continuous  time  quadrature  form 
by  the  equation 


c(t)  =  Re  jc(t)  ej27r£ot}  ,  (3) 

where  c(t)  is  the  continuous  time  form  of  Eq.  2  and  f  is  the 
carrier  frequency. 

Since  the  echo  signals  are  assumed  to  be  narrowband,  c(t)  is 

a  slowly  varying  function  with  respect  to  fQ.  If  the  target's 

location  is  restricted  to  the  main  lobe  of  the  transducer  beam 

patterns,  the  maximum  correlation  at  time  ^max  will  occur  within 

Tq/2  of  cCO),  and  cCt^)  *  c(.0)  . 

The  maximum  crosscorrelation  c(t  )  is  related  to  the  target 

max 

dimension  whereas  t  is  related  to  the  location.  This  will  occur 

max 

when  -2irfot  equals  the  angle  of  c(t)  in  Eq.  3.  Given  that  c(t)  is 
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slowly  varying,  |c(0)  |  =  cO^)  and  arg{c(0)}  =  -27^^. 
Accordingly,  the  sampled  version  of  c(0)  provides  sufficient  infor¬ 
mation  to  calculate  the  dimension  and  location  estimates. 

Expressions  for  the  horizontal  and  vertical  dimension  and 
angular  location  estimates  are  easily  obtained  for  the  case  where 
individual  fish  echoes  are  not  overlapping  (see  Appendix  A) .  The 

analysis  for  the  more  general  case  where  many  fish  contribute  to  a 

6 

given  quadrature  sample  is  given  by  Gardner  and  Jackson.  Both 
analyses  result  in  the  following  expressions  for  the  horizontal  and 
vertical  spatial  estimates. 


DH  ■  ^  V1-|£LRt0)l  f2''51  R 


Dv  =  y 


✓2 
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v/l-|CUD(0)|  (2/5)  R 


„  1,-1  ImfeLR(0)l 

6h"  -Y^ta" 


Re)cu(°)| 

■  1  ImfeuD(0)( 


%  =  '  YV  Re  (Cu0(0)(  ' 


(m) 

On) 

(deg) 

(deg) 


(4) 


where  R  is  the  target  range  and  yM  and  y,,  are  the  conversion  fac- 
tors  from  mechanical  to  electrical  degrees,  y^  =  (2irfoSpj)/c  and 
Yy  =  (2TrfQSy)/c,  where  is  the  separation  between  the  horizontal 
transducer  centers  (5.7  cm),  Sy  is  the  separation  between  the 
vertical  transducer  centers  (8.6  cm),  and  c  is  the  assumed  sound 
velocity  in  water  (1500  m/s) . 
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The  term  /2/y  /1-|2(0)|  represents  the  standard  deviation 
(with  respect  to  the  random  positions  of  the  scatterers)  of  the 
azmuthal  angle.  The  term  2/5  converts  the  standard  deviation  of  a 
collection  of  points  bounded  by  an  ellipsoid  to  the  dimension  of 
the  ellipsoid  (see  Appendix  B) . 

The  down-range  (x’  direction)  dimension  estimate,  D^,  is 
calculated  using  the  energies  in  tb  jadrature  samples  (see 
Appendix  C) .  The  down-range  dimension  is  given  by 

tt|Z(i)|2)2 

0  ~  1800  _i _  (s 

°DR  f  *  C5 

s  £|z(i)|4 

i 

where  f  is  the  sampling  frequency  and  the  summations  are  over  all 
quadrature  samples  of  a  transducer  signal. 
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CHAPTER  3 


DOCUMENTATION  OF  THE  SIMULATION  PROGRAMS 


In  this  chapter,  a  more  detailed  account  is  given  of  the 
system  of  programs  used  in  simulating  the  split-beam  system.  These 
programs  are  mostly  written  in  Extended  FORTRAN  language,  and  were 
used  at  the  facilities  of  the  University  of  Washington  Academic 
Computer  Center.  A  listing  of  the  programs  with  sample  inputs  and 
outputs  is  included  in  Appendix  D. 

The  computer  simulation  of  the  split-beam  method  is  accom¬ 
plished  using  the  system  of  programs  and  files  shown  in  Fig.  4. 

The  main  component  is  the  program  called  SNGLPNG  which  simulates  a 
single  ping  on  an  elliptical  fish  school  and  estimates  the  school's 
size  and  location.  SNGLPNG  obtains  data  needed  for  each  run  from 
two  sources.  The  input  file  contains  the  desired  parameters  for 
the  run,  including  information  on  the  fish  school  and  echo  sampling 
procedure.  The  results  of  previous  runs  which  can  be  used  in  the 
current  run  are  contained  in  a  second  input  file  from  tape  storage. 
If  no  previous  runs  have  been  made  on  a  particular  school,  then  the 
initial  inputs  to  SNGLPNG,  created  by  EFILL,  are  the  scatterer 
locations  within  the  ellipsoid.  The  results  of  the  simulation  for 
several  pings  are  than  analyzed  and  plotted  using  other  FORTRAN 
programs  with  Numerical  Plotting  System  routines. 
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Figure  4.  Split -beam  simulation  structure 

3.1  Description  of  EFILL 

The  starting  point  of- the  simulation  is  the  program  EFILL. 

This  program  simply  generates  random  locations  within  a  specified 
ellipsoid  and  stores  these  locations  on  tape.  Three  random  numbers 
corresponding  to  the  x,  y,  and  z  coordinates  of  each  fish  are  gen¬ 
erated  using  the  intrinsic  computer  function,  RANF .  This  function 
produces  independent,  uniformly  distributed  random  numbers  within 
the  range  0,1.  The  three  values  from  RANF  are  shifted  to  the  range 
(-1,1)  by  subtracting  0.5  and  doubling  each  value.  The  correspond¬ 
ing  location  is  tested  to  determine  whether  it  lies  within  a  unit 
sphere.  If  the  point  does  not  lie  in  the  sphere,  three  new  random 
numbers  are  generated;  otherwise,  a  fish  location  is  obtained  by 
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multiplying  the  x,  y,  and  z  coordinates  by  the  values  of  a,  b,  and 
c  used  in  the  following  expression  for  defining  the  ellipsoid: 


This  process  is  repeated  until  the  desired  number  of  fish  locations 
is  obtained. 

The  format  used  to  store  these  locations  on  tape  is  compatible 
with  future  outputs  of  SNGLPNG  which  contain  additional  information 
about  the  fish  such  as  scattering  amplitudes  and  Doppler  shifts. 
Since  EFILL  doesn't  calculate  this  information,  that  space  is 
filled  with  zeroes. 

3.2  Description  of  SNGLPNG 

The  program  SNGLPNG  is  composed  of  five  sections  as  shown  in 
Fig.  5.  The  information  gathered  in  the  input  section  determines 
the  program  flow  through  the  remaining  sections.  For  example,  if 
only  the  sampling  frequency  is  changed  from  a  previous  run,  the 
fish  school  data  are  taken  from  the  previous  run,  and  the  program 
skips  directly  to  the  echo  processing  section. 

Input  Section 

The  input  section  reads  data  from  the  input  file  and  tape 
storage.  The  desired  parameters  for  the  run  are  provided  by  the 
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Figure  5.  Block  diagram  of  SNGLPNG 
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input  file  and  a  summary  oi  the  parameters  used  in  a  previous  run 
is  found  on  tape  storage.  If  the  fish  school  is  the  same  in  both 
files,  the  program  skips  to  the  echo  processing  section. 

Fish  School  Section 

The  fish  school  section  performs  calculations  involving  the 
individual  scatterers.  Round-trip  propagation  delays,  a  random 
scattering  amplitude,  corrections  to  this  amplitude  to  account  for 
beam  patterns,  and  a  random  Doppler  shift  are  computed  for  each 
scatterer.  If  the  desired  location  determined  by  the  range,  R,  and 
bearing,  or  the  desired  orientation  determined  by  the  aspect 

angle,  0,  differ  from  those  of  the  previous  run,  then  new  propa¬ 
gation  delays  are  needed.  The  path  lengths  associated  with  these 
delays  are  determined  by  the  geometry  shown  in  Fig.  6.  The  total 
acoustic  path  lengths  for  the  left,  right,  up,  and  down  transducer 
centers  are  path  1  plus  path  5,  path  4,  path  3,  or  path  2,  respec¬ 
tively.  The  propagation  delays  are  these  lengths  divided  by  the 
assumed  speed  of  sound  in  water  (1500  m./sec) . 

New  amplitude  corrections  reflecting  each  scatterer' s  new 
position  in  the  beam  patterns  must  also  be  computed  if  the  school 
location  or  orientation  is  changed.  Expressions  of  the  form 

sin6  m9+C 

— q —  are  fitted  to  transducer  calibration  curves  from  an 

actual  split-beam  transducer.  The  results  of  these  fits  are  given 
in  Appendix  E. 


SIDE  VIEW 


z 


PATH  3  =  j[R-  X(I)]2  +  Y(I)2  +  [Z(I)  -  SV/2]2J 


TOP  VIEW 


PATH  5  =  j[R  -  X(I)  -  SHS]2  +  [Y(I)  +  SHC]2  +  Z(I)2J  ^ 
Positive  values  of  ip  are  to  the  left  of  the  beam  axis 
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Figure  6.  Illustration  of  propagation  path  lengths. 
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The  horizontal  and  vertical  angles  describing  the  angular 
location  of  each  fish  are  given  by 

THH  =  ip  -tan"1  [y/(R-x)]  (rad) 

and  (6) 

THV  =  tan-1  [  z//(R-x)2  +  y2  ]  (rad). 


where  x,  y,  and  z  are  the  fish  coordinates,  and  R  and  <p  are  the 
range  and  bearing  to  the  center  of  the  school. 

This  two-dimensional  location  is  converted  to  an  equivalent 
horizontal  or  vertical  angle  which  has  the  same  beam  strength  by 
the  equation 


THHEQ  =  /tTH2  +  (e*THV)2 

ir/0.439 

(rad) 

THVEQ  =  /(THH/e)2  +  THV2 

ir/0.291 

(rad) , 

where  e  is  the  horizontal  to  vertical  eccentricity  of  each  beam 
pattern.  THHEQ  is  calculated  for  the  transmit  and  up-down  beam 
patterns,  and  THVEQ  is  calculated  for  the  transmit  and  left-right 
beam  patterns.  The  factors  ir/0.439  and  tt/0.291  are  needed  to  make 
the  first  null  of  sin0/0  correspond  to  the  main  lobe  boundaries  of 
the  calibration  curves. 

THHEQ  or  THVEQ  can  now  be  used  in  the  equations  chosen  to  fit 
the  horizontal  or  vertical  calibration  curves  for  the  transmit, 
left-right,  and  up-down  transducer  sections.  These  equations 
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determine  the  beam  losses  of  the  transmit  signal  combined  with 
either  the  left-right  or  the  up-down  echo. 


where 


El 

=  0.017189  • 

E4  J 


THVEQS 

THVEQ 

THHEQS 

THHEQ 


THHEQS  and  THVEQS  are  computed  for  the  transmit  beam,  THVEQ  for  the 
left-right  beam,  and  THHEQ  for  the  up-down  beam. 

If  there  are  no  changes  in  the  school  location  or  orientation, 
then  the  propagation  delays  and  beam  corrections  from  the  previous 
run  are  read  from  tape  storage.  The  program  then  proceeds  to  the 
scattering  amplitude  calculations.  If  new  amplitudes  are  desired, 
they  are  generated  from  a  Rayleigh  process.  The  distribution 
function  for  a  Rayleigh  process  is  given  by  , 

y  =  P(x)  =1  -  e 

If  values  of  y  are  chosen  independently  from  a  uniform  distribution 
ranging  from  zero  to  one,  the  inverse  of  P  is 
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x  =  P_1(y)  =  ^2^~in~y  , 

and  the  corresponding  values  of  x  will  be  Rayleigh  distributed. 

The  computer  function,  RANF,  generates  such  values  of  y  for  each 
fish.  The  value  a,  which  is  proportional  to  the  mean  of  the 
Rayleigh  distribution,  is  also  generated  by  RANF  for  each  fish, 
and  is  shifted  to  range  uniformly  from  0.5  to  1.5.  The  result 
is  that  the  amplitude  for  each  fish  is  chosen  from  a  Rayleigh 
distribution:  however,  the  distribution  mean  is  chosen  indepen¬ 
dently  for  each  fish,  simulating  a  different  size  for  each 
fish. 

If  no  changes  were  needed  in  the  amplitudes,  the  previous 
amplitudes  are  read  from  tape  storage,  and  the  program  proceeds  to 
the  Doppler  shift  calculations.  The  Doppler  shifts  are  generated 
independently  from  a  uniform  distribution  ranging  from  20  to  100  Hz 
when  a  30  kHz  carrier  frequency  is  used.  Random  numbers  generated 
by  RANF  are  shifted  from  the  range  (0,1)  to  obtain  those  Doppler 
shift  values.  Note  that  the  Doppler  shifts,  amplitudes,  and  loca¬ 
tions  are  generated  independently  and  thus  are  uncorrelated. 

Echo  Processing  Section 

The  next  part  of  SNGLPNG  is  the  echo  processing  section.  The 
main  functions  of  this  section  are  shown  in  Fig.  3.  The  first  task 
is  the  formation  of  the  quadrature  samples.  The  propagation  delays 
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axe  scanned  and  the  times  of  the  first  and  last  samples  are  de¬ 
termined,  The  first  sample  is  computed  one  sampling  period  after 
the  arrival  of  the  first  echo,  and  the  last  sample  is  computed  at 
the  last  sampling  interval  prior  to  the  end  of  the  echo.  Signal 
and  noise  samples  are  calculated  at  the  desired  sampling  rate  until 
the  end  of  the  last  fish  echo.  If  the  signal  is  filtered,  addi¬ 
tional  samples  are  calculated  from  the  decaying  filter  output. 

In  order  to  achieve  the  desired  signal-to-noise  ratio,  the 
gain  or  magnitude  of  the  noise  samples  is  varied.  This  gain  factor 
is  determined  by  the  expected  values  of  the  signal  and  noise  ener¬ 
gies,  so  the  actual  SNR's  obtained  in  the  simulation  should  fluc¬ 
tuate  around  the  desired  value.  The  expected  signal  and  noise 
energies  are  initially  computed  for  the  case  where  no  filter  is 
applied.  For  a  rectangular  transmit  pulse,  the  expected  signal 
energy  is  the  sum  of  the  received  energies  (corrected  amplitude 
squared)  from  all  fish  times  the  number  of  samples  taken  of  each 
fish  echo.  It  is  assumed  that  the  incidences  of  constructive 
interference  between  fish  echoes  will  be  counterbalanced  by  inci¬ 
dences  of  destructive  interference. 

The  expected  energy  of  the  left-right  echoes  will  differ  from 
that  of  the  up-down  echoes  because  of  the  difference  in  beam  pat¬ 
terns.  The  expected  signal  energy  used  to  calculate  the  noise 
gain  factor  is  the  average  of  the  left-right  and  up-down  expected 
energies. 
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For  a  Hanning  transmit  pulse,  the  expected  energy  is  reduced 
from  the  rectangular  case  by  the  ratio  of  the  energy  in  the  two 
pulses.  The  energy  in  a  unit  amplitude  rectangular  pulse  is  equal 
to  the  pulse  width,  PW,  times  the  amplitude  squared. 

The  energy  in  a  Hanning  pulse  with  unit  peak  amplitude  is 
given  by 

E  =  —  f  sin4  t  dt  =  f  PW. 

IT  J  o 

O 

Therefore  the  expected  signal  energy  using  a  Hanning  pulse  is  0.375 
times  the  expected  energy  using  a  rectangular  pulse  of  the  same 
length. 

To  calculate  the  expected  energy  of  the  noise  samples,  the 
process  from  which  these  samples  are  generated  must  be  known.  The 
real  and  imaginary  components  of  each  noise  sample  are  independent 
identically  distributed  normal  random  variables  having  zero  mean  and 
unit  variance  (i.i.d.  N(0,1)).  The  method  of  generating  these  sam¬ 
ples  will  be  discussed  later.  The  expected  energy  in  a  noise  sample 
is  equal  to  the  expected  value  of  the  magnitude  squared  which  is 
the  sum  of  the  component  variances  since  the  magnitudes  have  zero 
mean  values.  The  total  noise  energy  expected  is  then  two  times  the 
number  of  samples. 

Given  the  expected  signal  energy  and  a  desired  SNR,  a  gain 
factor  is  computed  by  the  subroutine  NGAIN.  The  noise  samples. 
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when  multiplied  by  this  gain  factor,  will  have  an  expected  energy 
necessary  to  obtain  the  desired  SNR. 

The  next  step  in  the  echo  processing  .section  is  to  compute  the 

unfiltered  signal  samples.  This  process  is  also  described  in 

Chapter  2.  For  each  sample  time,  t^,  the  magnitude  of  the  echo 

envelope  is  computed  for  each  fish  by  the  function  PULSE.  If  the 

echo  from  a  fish  is  present  at  the  transducer,  the  value  of  the 

2 

envelope  is  equal  to  one  for  the  rectangular  pulse  or  sin  (t^-T) 
for  the  Hanning  pulse,  where  T  is  the  propagation  delay  to  the 
transducer-half  center.  If  the  fish  has  not  yet  arrived  at  or  has 
already  passed  the  transducer,  the  value  returned  by  PULSE  is  zero. 
A  nonzero  value  from  PULSE  causes  that  fish  echo  to  be  included  in 
the  sample  described  by  Eq.  1.  Additional  signal  samples  are 
computed  at  the  desired  sampling  rate  for  the  duration  of  the 
school  echo. 

An  equal  number  of  noise  samples  are  now  calculated  using  the 
RANF  function.  As  stated  earlier,  the  real  and  imaginary  compo¬ 
nents  of  the  noise  samples  are  i.i.d.  N(0,1)  random  variables.  To 
obtain  these  components,  two  uniformly  distributed  outputs  from 
RANF  are  used  in  the  equations 


R  =  /-2  Zn  [RANP(0.)J 
0  =  2tt  RANF(1.) 


26 


u  =  R  cos6 
v  =  R  sin0. 

The  noise  sample  is  composed  of  the  components  u  and  v  multiplied 
by  the  gain  factor  computed  in  NGAIN. 

If  the  signal  and  noise  samples  are  to  be  filtered,  then  the 
program  proceeds  to  the  filter  calculations.  A  single-pole  digital 
filter  is  used  whose  transfer  function  is  given  by 


H(z)  = 


1-a 


-1  ' 


1-az 


This  corresponds  to  the  iterative  equation 


y  =  (l-a)x  +  a  y  ,  , 
Jn  v  '  n  'n-1  * 


(9) 


where  x  and  y  are  the  input  and  output  of  the  filter,  and  the  sub¬ 
scripts  n  and  n-1  denote  the  n*h  and  preceding  values  of  the 
sequences.  Given  a  3  dB  bandwidth  frequency,  f^w>  the  value  of  "a" 
can  be  determined.  The  steady-state  transfer  function  is  obtained 


by  replacing  z  by  e 


H(e^wTs) 


1-a 


1-ae 


-jwTs 


At  the  3  dB  frequency,  the  magnitude  of  the  steady- state 
transfer  function  squared  is  0.5  by  definition,  so  the  value  of  "a" 
for  a  given  f.  is  found  by  solving 


0.5  = 


1-a 


l_ae“j27rfbw/fs 


(1-a)' 


2a  cos(2Trf^w/fs)  +  az 


(10) 
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or 

a  =  -cos(2Trfbw/fs)  +  2  -  /[2  cos(2irfbw/fs) -4] 2  -  4  /  2  . 

Values  of  ''a’1  range  from  0.172  to  1.0,  corresponding  to  values  of 

fbw  £rom  V2  t0  °* 

The  first  use  of  the  filter  is  to  compute  the  expected  energy 
of  the  filtered  signal  samples.  Again  the  assumption  is  made  that 
the  constructive  and  destructive  interferences  between  individual 
fish  echoes  cancel  each  other,  and  the  expected  energy  is  the  sum 
of  the  filtered  individual  fish  echoes.  The  integration  of  the 
filtered  rectangular  and  Hanning  pulses  is  performed  numerically  by 
summing  the  results  of  Eq.  9.  The  sampling  interval  is  the  pulse 
width  divided  by  10  for  the  rectangular  pulse  and  by  50  for  the 
Hanning  pulse.  The  process  is  stopped  when  the  energy  of  a  sample 
is  less  than  1%  of  the  sum  of  the  energies  of  the  previous  samples. 

The  expected  filtered  noise  energy  is  found  analytically  by 
finding  the  ratio  of  the  filtered  to  unfiltered  noise  energies  in 
the  frequency  domain.  The  frequency  spectrum  of  the  filtered  noise 
is  the  product  of  the  filter  spectrum  times  the  constant  spectrum 
of  the  unfiltered  white  Gaussian  noise.  By  Parseval's  theorem,  the 
integral  of  the  square  of  the  spectrum  magnitude  is  equal  to  2tt 
times  the  integral  of  the  signal  energy  in  the  time  domain.  Thus 
the  ratio  of  the  filtered  to  unfiltered  noise  energies  can  be  found 
by  integrating  each  spectrum  over  a  range  equal  to  one  period  of 
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the  filtered  spectrum  and  computing  the  ratio  of  the  results, 
integral  of  the  filtered  spectrum  is 

-  .2,  (l-a)2  2tt*Ts*  (l-a) 

s  JQ  1  +  a2  -  2a  cos0  1+a 


The 


where  the  variable  change  wT  -*■  0  is  made. 

The  integral  of  the  unfiltered  spectrum  involves  a  constant 
spectrum  of  unit  height, 

2tt 

Ts  /  lde.2,Ts, 

o 

and  the  ratio  of  unfiltered  to  filtered  expected  noise  energies  is 
(l-a)/Cl+a). 

Before  the  signal  and  noise  samples  can  be  filtered,  an  ini¬ 
tial  value  for  the  filtered  noise  just  prior  to  the  first  signal 
sample  must  be  computed.  The  expected  value  of  this  noise  sample 
is 

E[yQ]  =  Efay^  +  (l-a)xQ] 

=  E {a[ay_2  +  (l-a)x  j]  +  (l-a)xQ} 

2 

=  (l-a)  E[xq  +  ax_j  +  a  x_2  +  ***] 

=  (l-a)  E[x]  (1  +  a  +  a2  +  •••) 

=  E[x], 


E[x]  =  E[xq]  =  E[x_1]  =  •  • 


where 
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The  variance  of  yQ  can  also  be  expressed  in  terms  of  the  input 
variable  x. 

E[y^]  =  E{(l-a)2  [xQ  +  ax_1  +  a2x_2  +  •••]  } 

n  ,2cr2  22  .42  , 

=  (1-a)  E[xq  +ax1+ax2+  ••• ] 

=  (1-a)2  E[x2]  (1  +  a2  +  a^  +  •  ••) 

=  ETx2! 

1+a  tLX  J  » 


where 


and 


E[x  x  1=0  except  for  u  =  v 
n-u  n-vJ  1 


E[x2]  =  Etx2^]  =  E[x2_2]  = 


The  variance  of  yn_^  is  given  by 


=  E[y2]  -  E [y  ] 2 
J  o 


1-a  cr  2,  r  ,2 
E[x  ]  -  E [x]  . 


The  ratio  (1-a)/ (1+a)  is  also  the  ratio  of  the  unfiltered  to 


filtered  sample  energies,  so  a  reasonable  value  for  yQ  is  a  single 
unfiltered  noise  sample  multiplied  by  the  square  root  of  this 


ratio . 
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The  components  of  the  filtered  signal  and  noise  samples  are 
each  computed  using  Eq.  9.  After  the  last  signal  sample  is  used, 
the  values  of  xn  in  Eq.  9  for  the  signal  components  are  zero  while 
the  filtered  signal  decays.  Additional  noise  samples  are  generated 
and  filtered  during  the  signal  decay,  and  the  filtering  process 
continues  until  the  SNR  for  a  filtered  sample  drops  below  the 
desired  echo  SNR  divided  by  10.  At  the  end  of  the  filtering  proc¬ 
ess  the  program  skips  to  the  sample  correlation  and  normalizing 
calculations. 

The  part  of  the  program  skipped  above  is  used  on  runs  where 
only  the  desired  SNR  is  changed.  This  part  is  reached  only  if  no 
changes  were  made  in  the  fish  school  section  or  the  sampling  param¬ 
eters.  The  signal  and  noise  samples  from  the  previous  run  are 
used,  and  a  new  noise  gain  is  computed  by  NGAIN  from  the  new  de¬ 
sired  SNR.  The  previous  noise  samples  are  multipled  by  the  new 
gain  factor  and  the  program  proceeds  to  the  echo  crosscorrelation 
calculations. 

If  the  desired  SNR  is  infinite,  the  crosscorrelations  of  the 
left-right  and  up-down  echoes  are  computed  using  only  the  unnor¬ 
malized  signal  samples  in  Eq.  2.  For  finite  SNR's,  the  cross¬ 
correlations  are  computed  using  unnormalized  samples  of  signal  plus 
noise.  Next,  the  signal  and  noise  samples  are  normalized  by  the 
square  root  of  their  respective  echo  energies.  The  echo  cross¬ 
correlations  are  normalized  by  the  square  root  of  the  product  of 
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the  energies  of  the  two  correlated  echoes.  Other  values  relating 
to  the  energy  in  each  sample  and  to  the  summed  energy  of  all  fish 
contributing  to  each  sample  are  also  normalized.  These  values  are 
used  in  estimating  the  down-range  dimension  of  the  school  as  de¬ 
scribed  in  Chapter  2. 

The  final  part  of  the  echo  processing  section  is  the  calcula¬ 
tion  of  the  fish  school  location  and  dimension  estimates.  The 
normalized  left-right  and  up-down  crosscorrelations  are  used  in 
computing  the  horizontal  and  vertical  angular  location  and  dimen¬ 
sion  estimates  as  described  by  Eq.  4.  The  down-range  dimension 
estimate  which  is  related  to  the  sample  energies  is  given  by  Eq.  5 
or  Eq.  21. 

Data  Storage  and  Print  Sections 

The  results  of  the  simulation  are  now  stored  on  tape.  Two 
tapes  are  used  in  order  to  minimize  information  retrieval  time 
during  further  data  processing.  On  the  first  tape,  the  complete 
data  needed  for  future  runs  are  stored,  while  on  the  second  tape 
only  the  final  spatial  estimates,  echo  crosscorrelations,  and  echo 
energies  are  stored  for  use  in  analyzing  and  plotting  the  results,. 

The  final  section  in  SNGLPNG  is  the  print  section.  During  the 
execution  of  the  previous  sections,  a  summary  of  the  input  file 

i 

data,  tape  storage  parameters,  and  messages  denoting  the  program 
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flow  are  printed,  but  the  results  of  the  simulation  are  printed  in 
this  section.  The  results  are  divided  into  four  parts.  The  print 
code  supplied  by  the  input  data  determines  which  parts  will  be 
printed.  The  first  is  a  listing  of  the  fish  coordinates,  their 
propagation  delays,  and  their  amplitude  and  Doppler  characteris¬ 
tics.  The  second  part  is  a  listing  of  the  normalized  quadrature 
samples,  the  number  of  fish  contributing  to  each  sample,  and  the 
normalized  energy  density  in  each  sample.  The  third  part  printed 
is  a  listing  of  the  unnormalized  signal  and  noise  samples.  The 
noise  samples  are  always  computed  for  future  runs  even  though  they 
aren't  used  when  the  desired  SNR  is  infinite.  The  last  part,  which 
is  always  printed,  contains  a  summary  of  the  input  data,  the  signal 

and  noise  energies,  the  echo  crosscorrelations,  and  the  dimension 

<  1  * 

estimates  computed  from  the  ellipsoid  dimensions,  the  fish  coordi¬ 
nates,  and  the  quadrature  samples. 

After  leaving  the  print  section,  if  results  are  d.esired  for 
another  SNR,  the  program  skips  back  to  "the  echo  processing  section 
and  computes  new  results  by  changing  only  the  gain  of  the  noise. 
Otherwise,  the  program  returns  to  the  input  section  to  begin  on  a 


new  school. 


CHAPTER  4 


ANALYSIS  AND  RESULTS  OF  THE  SIMULATIONS 


Two  general  areas  relating  to  the  performance  of  the  split- 
beam  technique  are  studied  using  the  simulations.  Section  4.1 
investigates  the  effect  of  adding  uncorrelated  noise  to  the  target 
echo.  Since  the  dimension  estimates  are  calculated  from  the  signal 
crosscorrelations,  biases  toward  larger  values  are  expected  be¬ 
cause  of  the  decorrelating  effect  of  the  noise.  Section  4.2 
attempts  to  determine  the  relationship  between  the  target,  transmit 
pulse,  and  sampling  parameters  in  estimating  target  characteristics. 
Fluctuations  in  these  estimates,  which  are  most  noticeable  when  the 
signals  are  undersampled,  are  traced  to  terms  in  the  crosscorrela¬ 
tion  calculations. 

Three  sizes  of  schools  at  a  range  of  500  m  are  simulated.  In 
the  large  school,  the  ellipsoid  dimensions  in  the  x',  y',  and  z 
directions  shown  in  Fig.  1,  are  set  equal  to  50  m,  75  m,  and  25  m. 

In  the  middle  size,  they  equal  25  m,  37.5  m,  and  12.5  m,  and  in  the 
small  school,  10  m,  15  m,  and  5  m. 

Two  shapes  of  transmit  pulse  envelopes  are  used.  A  1  ms 

rectangular  pulse  envelope  is  compared  with  a  1.385  ms  Hanning 
2 

(sin  )  pulse  envelope.  The  Hanning  pulse  is  chosen  because  its 
high-frequency  energy  spectrum  is  much  lower  when  compared  with  the 
rectangular  pulse. 

The  pulse  length  of  the  Hanning  envelope  is  determined  by 
equating  the  pulse  resolutions  defined  by  Eq.  11  for  both  shapes. 
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00  ^ 

/  C^(T)dT  co 

X  =  — - ,  where  cs(x)  =  J  s(t)s(t-x)dt  .  (11) 

2  Cg(0) 

For 

s(t)  J1  0<t<lms 
( 0  otherwise 

or 

{sin2  t~7oE  t  0  <  t  <  1 . 385  ms 
1.385 

0  otherwise 

the  value  of  the  pulse  resolution  is.  T  =  0.333  ms. 

The  sampling  frequency  is  5  kHz  in  most  of  the  simulations. 

The  bandwidth  of  the  unfiltered  noise  samples,  which  is  determined 
by  the  sampling  frequency,  is  also  equal  to  5  kHz  in  most  cases. 

The  signal-plus-noise  samples  are  then  digitally  filtered  using  a 
1  kHz  bandwidth. 

Two  sets  of  beam  patterns  are  used  in  the  simulations.  The 
first  set  is  described  in  Chapter  3  and  Appendix  E.  The  second  set 
is  also  described  in  Appendix  E,  and  its  use  is  denoted  by  an 
asterisk  in  the  simulation  results. 

4 . 1  Effects  of  Uncorrelated  Noise 

Some  insight  into  the  effect  of  adding  uncorrelated  noise  is 
provided  by  analyzing  the  crosscorrelation  of  two  signals  plus 
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noise  from  a  single  point  target.  The  normalized  echo  crosscorre¬ 
lation  is  given  by 


Ci;2(0) 


00  *  rs 

f  j sOO+njCt)]  [s*(t)  e"-,0+n2*(t)]dt 
{ci,l(0)*c2,2(0)}1/2 


(12) 


where  c^  ^  and  c2  2  are  t^le  unnormalized  autocorrelations  of  the 
signals  and  are  used  to  normalize  the  echo  crosscorrelation,  s(t) 
is  the  complex  return  of  the  point  echo,  0  is  the  bearing  to  the 
target,  and  h^(t)  and  n2(t)  are  uncorrelated  complex  noise  returns 
at  centers  1  and  2.  By  defining 


=  J  s(t)s*(t)dt 


nl,l  =  /  n1(t)H1*(t)dt 
-00 


E 22  =/  n2(t)n2*(t)dt 

>  -OO 


00 

Eg  2  =  /  s(t)n2*(t)dt 
?  -00 


Eg  i  =  /  s(t)n1*(t)dt 
-00 


00 

/  n1(t)n2*(t)dt , 
-00 


the  crosscorrelation  is  expressed  as 


cl,2(°)  - 


ES,S  e  ^  +  ES,2  +  ES,1  e  +  El,2 


(13) 


yl (ES,S  *  2  H  SS,ll +  El,l)(BS,S *  2Re{e36Es>2}*  E^2 )  ] 


The  noise  terms  n^(t)  and  n2(t)  are  assumed  to  be  uncorrelated  with 
the  signal  s(t)  and  each  other.  These  signals  are  band  limited 
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and,  if  their  time-bandwidth  product  is  much  greater  than  unity, 
the  magnitudes  of  Eg  ^  Eg  2,  anc*  ^1  2  t>e  much  less  than 

^S,S’  ^1,1'  0r  ^2,2' 

Equation  13  can  be  expanded  using  the  binomial  theorem  for 
both  high  and  low  signal-to-noise  ratios.  For  the  high  SNR  expan¬ 
sion,  the  terms  involving  uncorrelated  factors  are  neglected  and 
the  crosscorrelation  becomes 


Ei,2«»  =  e'j6|1-sk 

where  Eg  g»E^  j~E2  2  and  SNR~Eg  g/tCE^  j  +  E2  2^2^  ‘ 


(14) 


This  equation  shows  that  the  location  estimate,  determined 
from  0,  should  remain  unbiased  over  this  range,  but  that  the  dimen¬ 
sion  estimate,  obtained  from  | 2(0)|,  Is  a  function  of  SNR"'1'. 

This  illustrates  the  decorrelating  effect  of  the  noise. 

For  the  low  SNR  expansion,  the  terms  involving  uncorrelated 
factors  are  retained  to  illustrate  the  result  when  the  SNR-K). 


SNR  + 


:S,1 


3S,2 


%2 


VE1,1  E2,2  /  VEl,i  E2)2  Ve1,1  E2,2 


(15) 


where  E^  *  B2>2»ES <s  . 

Equation  15  shows  that  if  the  SNR  is  much  greater  than  the 
terms  involving  uncorrelated  factors,  which  are  small,  the  location 
estimate  will  remain  unbiased;  however,  the  dimension  estimate  is 
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severely  biased.  As  the  SNR-K),  the  remaining  terms,  which  all 
have  random  angles,  will  cause  the  location  and  dimension  estimates 
to  fluctuate  randomly  from  ping  to  ping. 

The  results  of  the  simulations  for  the  three  school  sizes 
agree  with  the  single  point  analysis,  as  is  illustrated  in  Figs.  7 
through  10.  All  three  schools  are  located  at  a  horizontal  bearing 
of  12.5°.  The  numbers  of  fish  in  each  school  were  chosen  so  that 
approximately  the  same  number  of  fish  contribute  to  a  given  sample 
for  all  school  sizes.  Figure  7  shows  the  mean  horizontal  bearing 
to  the  schools  versus  SNR  for  30  simulated  pings  on  30  indepen¬ 
dently  generated  schools,  for  each  school  size.  Figure  8  shows  the 
standard  deviations  of  these  estimates  versus  SNR  over  the  30 
pings.  These  figures  show  that  the  location  estimates  remain 
unbiased  and  have  fairly  low  standard  deviations  down  to  SNR's  of 
approximately  -5  dB. 

Figures  9  and  10  show  the  relationship  between  the  horizontal 
dimension  estimates  and  the  SNR  for  the  three  school  sizes.  Each 
dimension  estimate  becomes  biased  from  its  actual  value  as  SNR 
decreases,  and  each  asymptotically  approaches  the  slope  predicted 
by  Eq.  14  for  a  single  point  target. 

The  down-range  dimension  estimates  are  not  computed  from  the 
echo  crosscorrelations.  Instead,  they. are  derived  from  the  approx¬ 
imate  radial  density  function  using  the  normalized  energies  of  the 
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Figure  7.  Illustration  of  the  effect  of  uncorrelated 
noise  on  the  horizontal  bearing  estimate 
(actual  bearing  =  12.5°) . 


SIGNAL  TO  NOISE  RATIO  (OBJ 


Figure  8.  Illustration  of  the  stability  of  the  horizontal 
bearing  estimate  versus  uncorrelated  noise. 


SIGNAl  TO  NOISE  RATIO  (DB) 


Figure  9.  Illustration  of  biases  in  horizontal  dimension 
estimates  caused  by  uncorrelated  noise  (actual 
values  are  in  parentheses) . 


SIGNAL  TO  NOISE  RATIO  (DB) 


Figure  10.  Illustration  of  the  fluctuations  in  the 
horizontal  dimension  estimates  caused  by 
uncorrelated  noise. 
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echo  samples  as  described  in  Appendix  C.  The  normalization  reduces 
the  effect  of  additive  noise  on  the  dimension  estimates  and  the 
biases  seen  in  the  horizontal  and  vertical  estimates  aren't  ob¬ 
served  in  the  down-range  estimates.  Noise  should  change  the  shape 
of  the  density  function,  tending  to  reduce  its  peak  and  make  the 
density  more  uniform.  The  value  of  D  for  a  uniform  density  func¬ 
tion  PR(r)  =  1/2  for  -l<r<l  and  zero  otherwise  is 

D  =  j^1  (1/2) 2  drj"1  =  2  , 

which  is  the  actual  width  of  the  distribution.  Because  the  value 
of  D  for  the  density  function  corresponding  to  an  ellipsoid  is 
(5/6) *2,  an  increase  of  1/6  is  expected  in  the  simulations.  This 
increase  can  be  seen  in  Fig.  11. 


Figure  11.  Effect  of  noise  on  the  down-range 
dimension  estimate. 
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In  the  simulation,  the  beginning  and  end  of  the  school  is 
known;  however,  in  an  actual  system,  the  length  of  the  echo  is 
determined  by  thresholds.  Ping-to-ping  fluctuations  in  the  value 
of  will  increase  with  decreasing  SNR  as  the  echo  becomes 
harder  to  detect. 

4.2  Sources  of  Fluctuations  in  the  Spatial  Estimates 

The  intent  of  this  portion  of  the  study  is  to  identify  some  of 
the  parameters  that  affect  the  performance  of  the  split-beam  sys¬ 
tem.  Ping-to-ping  fluctuations  in  the  dimension  estimates  can  be 
traced  to  fluctuations  in  the  terms  of  Eqs.  1  and  2.  In  addition 
to  the  random  aspects  of  each  fish  echo,  the  terms  in  Eq.  1  are 
determined  by  the  school  size  and  density  and  the  transmit  pulse. 
The  terms  in  Eq.  2  are  determined  by  the  school's  down-range 
dimension  and  the  sampling  frequency.  For  certain  combinations 
of  these  parameters,  larger  fluctuations  in  the  dimension  esti¬ 
mates  are  noticed. 

To  better  understand  the  possible  sources  of  the  fluctuations, 
it  is  helpful  to  expand  one  term  of  Eq.  2  as  a  function  of  the 
individual  fish  echoes.  Initially  we  shall  assume  that  the  same 
fish  contribute  to  both  samples,  and  their  echo  amplitudes  are  also 
the  same  for  both  samples.  The  resulting  expansion  for  n  fish  is 
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2*  = 
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|aje^  + 

aj  ejAl  + 
+  (n^-n-2 


•  ••  +  ane^nj  ja Al^  +  ■••  +  a^e”^ ^n“^n) j 

...  +  a2  ejA'>  +  a.a  ej  +  a  a  c'j  (<t>l -4>n  "A1 ) 

n  In  In 

terms  involving  other  possible  pairs  of  fish  echoes) 


(16) 


where  "a"  represents  the  echo  amplitude,  A  represents  the  phase 
difference  between  the  right  and  left  echo  determined  by  the 
angular  position  of  the  fish,  and  (j)  represents  a  random  angle 
determined  by  the  left  signal  propagation  delay  to  each  random 
location.  The  sum  of  the  terms  in  the  form  a^e^  has  a  phase  angle 
bounded  by  the  angular  extent  of  the  school,  since  A  represents  the 
angular  position  of  each  fish.  The  next  cwo  terms  represent  phasors 
whose  angles  are  random.  More  random  phasors  result  from  the  other 
pairs  of  fish  echoes.  If  the  number  of  fish  in  the  samples  is 
large,  these  random  phasors  should  have  a  small  mean  expected 
value.  However,  when  only  a  few  random  phasors  are  added  to  a 
similar  number  of  phasors  within  the  school  boundaries,  the  angle 
of  the  product  Z^Z*  may  lie  outside  the  school  boundary.  This  will 
decrease  the  magnitude  of  the  left-right  crosscorrelation  and 
produce  an  abnormally  large  dimension  estimate. 

Now  we  consider  the  case  where  the  echo  amplitudes  are  not 
slowly  varying  with  time,  and  the  left  and  right  echo  amplitudes  of 
a  fish  at  the  time  of  the  sample  may  differ.  This  will  occur  for 
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samples  taken  near  the  edges  of  a  fish  echo.  If,  for  example,  the 

amplitude  of  the  right  echo  for  the  first  fish  is  zero,  then  the 

terms  and  a,a  e-^  ^  in  Eq.  16  are  zero.  For  small 

1  1  n 

numbers  of  fish,  the  loss  of  a  random  phasor  will  tend  to  increase 
the  fluctuations. 

Simulation  results  comparing  the  horizontal  dimension  esti¬ 
mates  for  schools  of  varying  densities  are  shown  in  Table  1.  The 
trend  toward  higher  standard  deviations  for  lower  fish  densities  is 
caused  by  the  occasional  occurrence  of  abnormally  large  estimates 
as  high  as  40  m  for  the  small  school.  Using  a  Hanning  pulse  whose 
amplitude  envelope  is  slowly  varying,  the  fluctuations  are  smaller 
than  for  the  rectangular  pulse,  although  the  trend  toward  increased 
fluctuations  for  decreasing  densities  is  still  noticed.  The  re¬ 
sults  of  increasing  the  sampling  rate  compared  with  changing  the 
pulse  shape  are  shown  in  Table  2.  If  the  number  of  samples  per 
pulse  length  is  increased  while  using  a  rectangular  pulse,  the 
fluctuations  caused  by  samples  taken  near  the  edge  of  the  fish  echo 
should  be  reduced  by  more  good  samples  taken  well  within  the  edges 
of  the  echo. 

The  above  discussion  suggests  that  the  maximum  bandwidth  of 
the  transmit  pulse  or  received  signals  prior  to  the  onset  of  unde¬ 
sirable  fluctuations  is  determined  by  the  spatial  characteristics 
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Table  1.  Illustration  of  dimension  estimate  fluctuations  versus 
fish  school  density  and  transmit  pulse  shape.  Results 
are  for  30  pings  on  independent  schools. 


Mean 

Horizontal 

Standard 

Simulation 

Number 

Dimension, 

Deviation, 

Parameters 

of  Fish 

Dh  (m) 

ODH  On) 

Small  school* 

500 

15.57 

3.58 

No  random  amplitudes 

No  Doppler  shifts 

250 

14.90 

2.50 

Rectangular  pulse 

100 

15.36 

5.66 

50 

13.30 

2.10 

25 

16.88 

8.99 

10 

14.61 

7.53 

Small  school* 

100 

14.28 

2.21 

Random  amplitudes 

Random  Doppler  shifts 

50 

13.69 

2.76 

Hanning  pulse 

25 

13.60 

3.62 

10 

12.21 

4.54 

Large  school*  500 

No  random  amplitudes  __n 

No  Doppler  shifts 

73.14 

72.54 

5.78 

5.23** 

Rectangular  pulse  125 

74.27 

6.04 

50 

71.86 

13.01 

Sampling  frequency  =  5  kHz 

School  locations:  ip  =  0,  range  =  500  m 

SNR  =  infinity 


*Alternate  beam  amplitude  corrections  used 

**Statistics  based  on  17  pings 
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Table  2.  Illustration  of  the  reduction  in  relative  dimension 
estimate  fluctuations  for  larger  school  sizes  (x* 
dimension)  versus  transmit  pulse  length.  Results 
are  for  30  pings  on  independent  schools  for  each 
school  size. 


Transmit 

Pulse 

Mean 

Horizontal 
Dimension 
Estimate, 
Dh  (m) 

Standard 

Deviation, 

0dh  M 

Sampling 

Frequency 

(kHz) 

Rectangular* 

15.36 

5.66 

5 

1  ms 

14.78 

3.27 

25 

Hanning* 

(sin2) 
1.385  ms 

13.78 

2.41 

5 

School  size:  small  (10  m,  15  m,  5m),  100  fish 
School  location:  \p=0°,  range  500  m 
No  random  amplitudes  and  Doppler  shifts  assigned 
to  scatterers 
SNR  =  infinity 

Alternate  beam  amplitude  corrections  used 


of  the  transducer.  The  rise  time  of  the  echo  amplitude  must  be  long 
compared  to  the  maximum  difference  in  the  arrival  times  of  a  fish  echo 
at  the  left  and  right  halves.  If  we  choose  the  rise  time  to  be  ten 
times  the  maximum  difference,  the  relationship  is  given  by 

At  SH  .  fl 
10  ~  1500  sinwmax  > 

where  At  is  the  100%  rise  time  of  the  echo,  SH  is  the  horizontal 
separation  of  the  transducer  halves  shown  in  Fig.  6,  and  6max  is 
the  maximum  practical  angle  in  the  main  lobe  for  receiving  a  fish 
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echo.  If  the  value  of  0  =  15°,  then  At  =  0.098  ms  and  the 

max  * 

bandwidth  of  the  signal  should  be  less  than  1/At,  or  10  kHz. 

Using  this  bandwidth,  the  sampling  rate  is  determined  only  by  the 
desired  number  of  samples  in  the  total  echo. 

Equation  16  can  also  be  used  to  illustrate  the  effect  of 
including  random  amplitudes  and  Doppler  shifts  on  the  dimension 
estimate  fluctuations.  The  random  Doppler  shifts  can  be  included 
in  A  by  redefining  A  as 

A  =  ¥  (fo+£d}  Sin0  » 

where  f^  is  the  Doppler  shift.  If  f0  =  30  kHz  and  the  largest 
value  of  f^  is  100  Hz,  then  the  change  in  A  caused  by  fj  is 
negligible  and  should  not  affect  the  dimension  estimates. 

The  effect  of  the  random  amplitudes  is  illustrated  by  com¬ 
paring  the  expected  amplitudes  of  the  random  phasors  for  the  equal 
amplitude  and  random  amplitude  cases.  When  the  amplitudes  are 
Rayleigh  distributed  random  variables,  the  parameter  a  of  the  dis¬ 
tribution  can  be  determined  by  normalizing  the  expected  echo  ener¬ 
gies  in  the  two  cases.  For  the  equal  amplitude  case,  using  a 
rectangular  pulse  and  neglecting  beam  amplitude  corrections,  the 
expected  echo  energy  is  given  by 

N  o 

k  l  AT  =  k  N  , 
i=l  1 


where  k  is  the  number  of  samples  per  pulse  length  and  A^  =  1  for 
all  fish.  For  the  Rayleigh  distributed  amplitudes,  the  expected 
echo  energy  is 


N  ? 
k  Z  AT 
i=l  1 


k  NE 


A2 

1 


kN  2a2. 


Equating  the  two  expected  energies  results  in  a  =  1//2. 

The  ratio  of  the  normalized  amplitudes  of  the  random  phasors 
is  given  by 


E[aiaj] 

E[l*l] 


IT 


4  ' 


The  fluctuation  in  the  crosscorrelation  caused  by  these  random 
phasors  should  be  reduced  for  the  random  amplitude  case  by  approx¬ 
imately  this  value. 

Since  the  dimension  estimates  are  a  function  of  the  square 
root  of  the  crosscorrelation,  their  fluctuations  should  be  reduced 
by  approximately  /jt/4.  This  reduction  should  be  statistically 
insignificant  for  the  simulations  shown  in  Table  3  if  the  esti¬ 
mates  are  assumed  to  be  normally  distributed. 

If  it  is  assumed  that  the  30  dimension  estimates  are  normally 
distributed,  an  approximate  value  for  the  variance  in  the  dimen¬ 
sion  fluctuations  can  be  computed.  The  expected  sample  variance 


is  defined  as 
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v  is  a  chi-square  statistic  with  29  degrees  of  freedom,*  and 

2 

the  corresponding  fluctuation  in  the  sample  variance,  a—  ,  is 

(2*29)/ (302)  O4  . 

_  2 

A  typical  value  of  v  for  the  small  school  is  9  m  ,  and  the  corre- 

2 

sponding  value  of  a—  is  2.4  m  .  These  fluctuations  are  greater 
than  the  reduction  of  /ir/4  caused  by  random  amplitudes. 


Table  3.  Illustration  of  the  effect  of  ascribing  random  ampli¬ 
tudes  and  Doppler  shifts  to  the  point  scatterers. 


School 

Size 


Small* 

(25  fish) 

Small* 

(100  fish) 

Medium 
(250  fish) 


Equal  Amplitudes 
No  Doppler  Shifts 


Mean 

Horizontal 

Standard 

Dimension, 

Deviation 

Dm  (m) 

Ohh  (m) 

15.38 

3.35 

13.79 

2.41 

36.33 

4.27 

Random  Amplitudes 
and  Doppler  Shifts 


Mean 

Horizontal 

Dimension, 

Dh  (m) 

Standard 
Deviation 
<7DH  (m) 

13.60 

3.62 

14 . 28 

2.21 

37.72 

4.61 

Transmit  pulse:  1.385  ms  Hanning 
Sampling  frequency:  5  kHz 
School  location:  \p  =  0,  range  =  500  m 
SNR  =  infinity 

*Alternate  beam  amplitude  corrections  used 
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Another  effect  on  the  dimension  estimates  is  the  expected 
reduction  in  fluctuations  as  the  pulse  length  is  shortened,  be¬ 
cause  higher  pulse  resolution  provides  more  statistically  inde- 

2 

pendent  echo  samples.  However,  if  the  pulse  length  is  shortened 
to  the  point  where  only  a  few  fish. contribute  to  each  sample,  then 
fluctuations  will  increase  as  described  earlier.  Simulation 
results  for  pulse  lengths  where  the  typical  number  of  fish  per 
sample  is  greater  than  10  are  shown  in  Table  4.  These  results 
illustrate  the  reduction  in  relative  fluctuations  for  lower  ratios 
of  pulse  resolution  to  down-range  size  of  the  school.  The  rela¬ 
tive  fluctuations  are  the  standard  deviations  divided  by  the 
actual  horizontal  dimension. 

The  distinction  between  single  and  multiple-fish  targets  can 
be  accomplished  using  the  ping-to-ping  fluctuations  in  the  cross¬ 
correlations.  If  two  fish  are  used  in  Eq.  16,  then  two  phasors 
pointing  in  the  direction  of  the  fish  will  be  added  to  two  random 
phasors.  If  the  separation  between  the  fish,  denoted  as  A^-A2»  is 
small,  then  the  two  random  phasors  of  equal  magnitude  will  be  in 
nearly  opposite  directions  and  very  little  fluctuation  will  be 
noticed.  However,  as  the  separation  increases,  the  random  phasors 
will  not  cancel  and  fluctuations  will  result.  A  single  fish  will 
result  in  no  fluctuations,  since  no  random  phasors  are  involved. 
Thus  the  detection  of  a  single  fish  at  a  given  range  can  only  be 
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accomplished  in  situations  where  the  fish  are  separated  by  an 
amount  sufficient  to  generate  detectable  fluctuations. 


Table  4.  Illustration  of  high  standard  deviations  of  dimension 
estimates  for  a  rectangular  transmit  pulse  and  low 
sampling  frequency. 


School 

Size 

Mean 

Horizontal 
Dimension, 
Dh  (m) 

Relative 

Standard 

Deviation, 

dh  h 

School  Size 
in 

x'  Dimension 
(m) 

Small* 

14.28 

0.155 

10 

Medium 

37.72 

0.122 

25 

Large* 

72.20 

0.090 

50 

2 

Transmit  pulse  envelope:  1.385  ms  Hanning  (sin  ) 

Sampling  frequency:  5  kHz 
School  location:  *p  =  0,  range  =  500  m 
Includes  random  amplitudes  and  Doppler  shifts 
for  scatterers 
SNR  =  infinity 

Number  of  fish:  small,  100;  medium,  250;  large,  500 
* 

Alternate  beam  amplitude  corrections  used 
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CHAPTER  5 


CONCLUSIONS 


Computer  simulations  of  the  split-beam  method  show  that  this 
technique  should  be  useful  in  estimating  the  spatial  characteris¬ 
tics  of  fish  schools.  The  dimension  estimates,  however,  are 
biased  by  noise  such  that  for  low  signal-to-noise  ratios  the 
dimensions  are  irretrievable.  For  example,  with  a  30  dB  signal- 
to-noise  ratio,  schools  as  small  as  15  m  can  be  estimated  at  a 
range  of  500  m;  however,  with  a  10  dB  ratio,  the  minimum  estimate 
obtainable  is  approximately  80  m.  Thus  all  schools  with  dimen¬ 
sions  of  around  80  m  or  less  would  have  the  same  dimension  esti¬ 
mate  at  this  signal-to-noise  ratio.  The  location  estimates  are 
not  biased  by  noise,  and  useful  location  information  should  be 
obtainable  for  signal-to-noise  ratios  down  to  -5  dB. 

Ping-to-ping  fluctuations  in  the  spatial  estimates  are  rea¬ 
sonably  small  for  sampling  frequencies  above  the  Nyquist  rate  and 
for  high  signal-to-noise  ratios.  These  fluctuations  decrease 
relative  to  the  school  size,  as  the  ratio  of  pulse  resolution  to 
down-range  extent  of  the  school  decreases.  Typical  results  give 
standard  deviations  of  approximately  2.5  m  for  the  15  m  school  and 
6  m  for  the  75  m  school.  No  changes  in  the  magnitude  of  the 
fluctuations  are  noticed  when  random  amplitude  and  Doppler  shifts 
are  added  to  the  point  scatterers.  A  small  reduction  in  the 
fluctuations  is  expected  for  the  addition  of  random  amplitudes. 
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but  the  sample  size  of  30  pings  is  not  large  enough  for  the  change 
to  be  noticeable. 

A  beneficial  use  of  the  fluctuations  is  in  the  detection  of 
single  fish.  This  can  be  accomplished  if  the  fish  in  a  given 
range  window  are  separated  by  the  angular  distance  necessary  to 
produce  detectable  fluctuations.  The  situation  will  still  exist 
where  two  smaller  fish  that  are  close  together  will  be  indistin¬ 
guishable  from  a  larger  fish. 

Undersampling  can  create  larger  fluctuations  for  certain 
combinations  of  the  target  size  and  density,  the  bandwidth  of  the 
transmit  pulse,  and  the  sampling  characteristics.  These  fluctua¬ 
tions  are  caused  by  random  terms  in  the  crosscorrelation  calcula¬ 
tions,  and  can  be  minimized  by  reducing  the  magnitudes  of  these 
random  phasors  and  including  enough  of  them  that  the  effect  of 
their  sum  is  small.  This  involves  increasing  the  number  of  fish 
per  sample  and  decreasing  the  sample  time-bandwidth  product. 

This  thesis  is  a  basic  study  of  the  factors  affecting  split- 
beam  spatial  estimates.  It  should  provide  useful  information  for 
actual  split-beam  measurements  and  implementations  of  the  split- 
beam  technique  in  larger  simulation  programs.  Suggestions  for 
further  research  involve  analyzing  actual  split -beam  data  in  an 
attempt  to  isolate  the  effects  studied  here. 
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APPENDIX  A 


DERIVATION  OF  SPATIAL  ESTIMATES  FOR  NONOVERLAPPING  ECHOES 


In  this  appendix,  the  expressions  for  the  angular  location  and 
dimension  estimates  are  derived  for  the  case  where  a  rectangular 
transmit  pulse  is  used  and  the  individual  fish  echoes  do  not 
overlap.  In  this  case  either  one  or  no  fish  contribute  to  each 
quadrature  sample  summation  in  Eq.  1.  Correspondingly,  each  non¬ 
zero  term  in  the  summations  of  Eq.  2  can  be  expresssed  as 
|z^|^  e^i,  where  is  the  normalized  magnitude  of  the  single  fish 
echo  present  at  time  t^,  and  <(k  is  the  angular  position  of  that 
fish  in  electrical  degrees.  The  summations  in  Eq.  2  over  all 
samples  can  be  expressed  as  summations  over  all  fish  by  knowing 
the  number  of  samples  per  fish  echo. 


ELR  -  l  |Zi|2  ^  ■  k  l  |Z»|2  ^  <17> 

where  i  is  summed  over  all  quadrature  samples,  n  is  summed  over  all 
fish,  and  k  =  fs  x  the  transmit  pulse  length. 

Now,  we  define  the  mean  and  variance  of  the  distribution  of 
angles  that  are  weighted  by  the  fish  sample  energies  as 


<J>  = 


Z  1 
n  1 

M2 

and 

a2  - 

z 

n 

N 

2  <♦, 

Z 

n 

N2 

V 

Z 

n 

N 

T  N  2 


n 


(18) 


If  the  fish  are  concentrated  around  the  mean  angle,  (J>, 
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a  =  <j>  -  d>  «  1  and 

n  rn  * 


Z  I Z  I2  o2 
a2  =  n  I  nl  n 

(J)  v  1 7 :  |2 


Using  the  Taylor  series  expansion,  the  crosscorrelation  becomes 

L  =ke^E  I  Z  1 2  e^11 
LR  n  I  n| 

2 


ZkeHl  |Znf  »j*n  -T>  * 


The  values  of  from  Eq.  17  are  normalized,  so  by  definition 


?  |z.  2  =  k  z  zl2  =  1. 

i  I  1  n  nl 


From  Eq.  19,  we  have 
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n 

Z  1 2  cr 
n  j  n 

.  5  2» 2  ♦» 
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n 
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z„  2 

z  z  2 
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n  n 
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n 

Thus  Z  I Z  1 2  a  =0,  and  the  crosscorrelation  is 

nhl  " 

r  —  f I  A  n 2  \ 


q  =  n  -  —  o2  ) 

lLR  u  ■  2  V  ’ 


The  conversion  from  electrical  to  mechanical  degrees  is  accom¬ 
plished  by  using  the  transducer's  spatial  characteristics.  Figure  12 
shows  the  transducer  halves  separated  by  SH  and  a  target  located  by 
r  and  0.  If  r  »  SH  and  0  is  small,  the  difference  in  path  lengths 
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to  each  half  is  SH  sin0  ~  SH  0.  The  electrical  difference  in 
o>o  too 

radians  is  <j>  *  —  SH  0  =  y^9>  where  y^  =  —  SH  is  the  conversion 
from  mechanical  to  electrical  degrees. 


SH 


SH  sin# 

Figure  12.  Conversion  from  mechanical  to  electrical  degrees. 


In  mechanical  degrees,  the  crosscorrelation  is 
jYH^ 


CLR  ~  6 


1  2  2. 

(i  -  2  YH  O0)  . 


The  expression  for  the  mean  horizontal  angular  position  of 
the  fish  is 


0„  =  —  tan"1 
H  YH 


Im  C 


LR 


Re  C 


LR 


and  the  horizontal  dimension  of  the  school  is  given  by 


dh  -  •  R  •  F  -  f  f7 


'LR 


R  •  F 


where  F  is  the  conversion  from  the  standard  deviation  of  a  distri¬ 


bution  to  the  actual  dimension  of  the  distribution. 


The  vertical  angular  location  and  dimension  estimates  are 


given  by 


and  SV  is  the  separation  between  the  vertical  transducer  halves. 


APPENDIX  B 


RELATIONSHIPS  BETWEEN  THE  SCHOOL  DIMENSIONS  AND  THE 
•  STANDARD  DEVIATION  OF  THE  FISH  LOCATIONS 


In  this  appendix,  the  factor  is  calculated  to  convert  the 
standard  deviation  of  a  distribution  of  scatterers  to  the  actual 
size  of  the  group.  The  simulation  uses  an  ellipsoid  filled  with 
scatterers  whose  locations  are  independent  and  uniformly  distri¬ 
buted.  The  standard  deviation  of  this  distribution  of  scatterers 
is  computed  here. 

Since  the  density  of  scatterers  is  uniform  throughout  the 
ellipsoid,  the  probability  density  function  (PDF)  of  the  scatter¬ 
ers  is  equal  to  the  normalized  cross-sectional  area.  For  an 
ellipsoid  defined  by 


2  2  2 
x  y  z  1 

2  +  ,2  +  2  =  1  > 

a  b  c 

the  cross-sectional  area  parallel  to  the  x-z  plane,  for  example, 
is  given  by 


A(y)  =  irac  [1  -  (y/b)  ] 


-b  <  y  S  b  . 


The  normalizing  factor  is  giyen  by 


r  A(y)  dy  =  4/3nabc 
d-h 


and  the  PDF  is  expressed  as 


p  (y)  =  irac [1- (y/b)2]  _  3(b2-y2) 
^  4/3  7rabc  4b^ 


i 


4/3  7r  abc 


60 


The  variance  in  the  y  dimension  is  given  by 


°Y  =  E(y2}  -  E{y}2  = /7  y2^fldy 
7  J-b  4  b3 


2  2, 


So  the  relationship  between  the  actual  dimension,  2b,  and  the 
standard  deviation,  o  ,  is 

’  y 


2b=2/5a=F  a  and  F=2^5. 

y  y 


APPENDIX  C 


DERIVATION  OF  THE  DOWN-RANGE  DIMENSION  ESTIMATE 

The  down-range  dimension  estimate,  DpR ,  is  calculated  using  the 
definitions 

1 

D  =  - — 

/CO  0 

PR(r)z  dr 

0 

and 

Ddr  =  6/5  D  ,  (20) 

where  PR(r)  is  the  radial  density  function  of  the  fish.  If  the 

density  function  for  a  uniformly  filled  ellipsoid,  described  by 
2  2  2  2  2 

x  +y/b  +z/c  =1,  is  used  to  compute  DnD,  the  value  differs  . 

UK 

from  the  actual  down-range  width  of  the  ellipsoid.  For  this  case, 

n  cVl-x2  _  3(l-x2) 

Xw  4/3iibc  '  "  4 

D  *  ■  (n?  l  a-x2)2^  1  - 1  ■ 

Since  the  actual  dimension  equals  2,  the  correction  factor  of  6/5  in 
Eq.  20  is  needed  for  the  ellipsoids  used  in  the  simulations. 

Two  methods  of  approximating  the  radial  density  function  are 
used  in  the  simulation.  Both  involve  energies  associated  with  each 
quadrature  sample.  The  first  method  forms  the  density  function 
from  the  normalized  energy  of  the  samples. 
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The  result  in  terms  of  samples  is 


•  1 2 


P(i)  a  — 


E  I Z  Ci) 1 2  A 


and  the  value  of  the  down-range  dimension  is 


DDR  2  6/5 


izci)  i2y 


1800  'i 


(f  W)‘ 


;m4 


S  |Z(i)| 'A  rs  E  |Z(i)| 
i  i 


where 


A  =  =  (m/smpl)  . 


The  second  method  forms  the  density  function  by  neglecting  the 
phase  of  the  fish  echoes  and  summing  the  energy  magnitudes  of  all 
contributing  fish  echoes.  These  sums  are  also  normalized,  and  the 
density  function  is  given  by 

l  (A.(k)  B(k)j 2 

P(i)  =  — - - o—  , 

S  (E  A.(k)  B(k)\2  A 

whore  the  notation  of  Gq.  1  is  followed.  The  down-range  dimension 
estimate  is  then  given  by 
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APPENDIX  D 


LISTINGS  OF  COMPUTER  PROGRAMS 


This  appendix  provides  listings  of  the  programs  used  in  the 
simulations.  The  first  listing  is  the  program  EFILL  which  gen¬ 
erates  the  random  locations  of  the  point  scatterers.  The  second 
listing  is  of  the  main  program  SNGLPNG. 
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Table  5.  Listing  of  I3PILL 


c 

c 

c 

c 

c 


PROGRAM  EFILLI INPUT, OUTPUT, TAPE1 , TAPE5= I NPUT , TAPE6=OUTPUT) 
THIS  PROGRAM  CREATES  AN  INITIAL  FILE  FOR  THE  PROGRAM  SNQLPNG. 
THE  FILE  CONTAINS  N  FISH  LOCATIONS  INSIDE  AN  ELLIPSOID  WITH 
X,Y,  AND  Z  DIMENSIONS  OF  2*A,2*B,  AND  2*C  METERS.  SPACE  FOR 
OTHER  DATA  USED  BY  SNQLPNG  BUT  NOT  COMPUTED  HERE,  IS  FILLED 
WITH  ZEROES. 

READ C 5 , 500)  N.A.B.C 
500  FORMAT! 15 ,5X,3F10.0) 

DATA  (NF1  1 1)  ,  1=2,3)  ,  (RMS(I)  ,  1=1  ,500)  ,  (DUM6M) ,  I  si  ,6)  , 

#  (DUM29I I) , 1=1 ,28) , (DUM8II) , 1=1 ,B)/-1 ,-1 ,500*1 . ,43*0/ 
DATA  (DUM500I  I)  ,  1=1  ,500/500*0/ 

J  =  1 


READ ( 5 ,501 )  NF1I1) 

CALL  RANSET (NF1  (1) ) 

DO  100  1=1,200 
100  DUMP=RANF  (1  . ) 

GOTO  20 

10  READ (5 , 501 )  NF1  (1) 

501  FORMAT! 13) 

I F (EOF (5) . NE . 0 . )  999,20 
20  WRITE (1)  NF 1 , DUM6 , DUM6 , DUMB , DUM6 
WRITE (1)  DUM29 

WRITE (6, 600)  (NF1 II) , 1=1 ,3) ,N,A,B,C 
600  FORMAT!///, IX, *RANDOM  FILLING  OF  ELLIPSOID  AND  TAPE  GENERATION* 

*  ///3X,*SCH00L  NUMBER  =* , 1 3 , 2  1 2 , 1 0X , 

*  //3X,*NUMBER  OF  FISH  =*,I6,/3X, 

*  *A IX  AXIS)  =* , FI 0 . 2 , 5X ,*B ( Y  AXIS)  =*,F10.2,5X,*CIZ  AXIS)  =*, 

*  FI  0 .2 , //5X , MNO.M , 6X ,*X*,9X ,*Y* , 9X , *Z* , //) 

DO  200  1=1 ,N 

90  0=2 .* (RANF  !0 . ) - . 5) 

P=2 .* (RANF (0 . ) - . 5) 

Q=2 .* (RANF (0. ) - . 5) 

TEST=0**2+P**2+Q**2 
IF ITEST.GT. 1 . )  GOTO  90 
X ( I ) =A*0 


Y ( I ) =B*P 
Z ( I ) =C*Q 

200  WRITEI6.610)  I ,X 1 1 ) , Y ( I ) ,Z ( I ) 

610  FORMAT (IX, 16.3F10.2) 

IF ( H . GE .500)  GOTO  190 
K=N+1 

DO  195  I =K , 500 
195  X ( I ) =Y 1 1 ) =Z ( I ) =0 . 

190  WRITE II  X , Y , Z , DUM500 , DUM500 , DUM500 , DUM500 , DUM500 , DUM500 , 
*  RMS , DUM500 , DUM500 , DUM500 , DUM500 
J=J  +  1 

IFU.GT.30)  GOTO  999 
GOTO  10 
999  STOP 
END 
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Table  6.  L.i sting  of  SNGLPNG 


PROGRAM  SNGLPNGI INPUT, OUTPUT, TAPE1 , TAPE2 , TAPE4 , TAPE6= INPUT, 

#  TAPE6=0UTPUT) 

C######################################################################* 


# 

THIS  PROGRAM  SIMULATES  A  SINGLE  PING  ON  A  FISH  SCHOOL  # 

AND  CALCULATES  THE  SCHOOL’S  LOCATION  AND  DIMENSION  # 

ESTIMATES  USING  THE  SPLIT-BEAM  METHOD.  THE  PROGRAM  # 

REQUIRES  TWO  INPUT  FILES:  # 

INPUT  FILE  -  DESIRED  PARAMETERS  FOR  THIS  RUN  # 

TAPE1  -  RESULTS  FROM  A  PREVIOUS  RUN  OR  FROM  # 

THE  PROGRAM  EFILL  H 

AND  PROVIDES  TWO  BINARY  OUTPUT  FILES:  * 

TAPE2  -  COMPLETE  RESULTS  FOR  USE  IN  FUTURE  RUNS  # 

TAPE4  -  SUMMARY  OF  RESULTS  FOR  DATA  ANALYSIS.  # 

# 

FURTHER  INFORMATION  IS  GIVEN  IN  SEPERATE  DOCUMENTATION.  # 

# 


#################################################*##################### 
DIMENSION  AMP (500) , ASNR (4) , A8NR1 (4) , BCL (500) .BCUI500) ,DHQ(2) , 

#  DRQ  f  2 ) .EFLI500) .EFU1500) ,ESI4) ,EN (4) ,EQ(4) . ESN (4) , ENK ( 4 ) , EQN (4) 

#  , ESS  1 500 , 4) ,RMS (500) , NF (3) , NP (3) , NS (4) ,NSEED(99) ,NF1  (3) ,NS1  (4) , 

#  NPWMOO)  ,0(100,4,2)  ,  QS  ( 1 00 , 4 , 2 ) ,  ON  MOO,  4, 2)  ,  SSNR  1100,4)  , 

#  SDQM00.4)  ,SDA 1 100,2)  ,SF MOO)  ,T (500,4)  ,W(4)  ,X (500)  ,X2(4)  , 

#  Y (500) , Z (500) , DVQI2) , DQ(4) ,SDAN(100,2) , DSN  10 ) 

COMMON  /PLSE/PW , FO , DP  1 500 ) /ANGL/A ,A2,B2,PI 

CONSTANTS  AND  INITIAL  VALUE 


IR=1 

Pl=3. 141592654 
TWOP I =2 .NP I 
CW=1 500 . 

RADEG=57. 2957795131 
SR2=SQRT (2 . ) 
SR5=2#SQRT(5.) 


MMM I NPUT  SECTION### 


TABLE  OF  99  RANDOM  SEEDS 

READ (5 , 5010)  (NSEEDI 11,1=1,99) 

5010  FORMATI8F10.8) 

NDESIRED  SNR  VALUES  FOR  ADDITIONAL  RUNS# 

READ (5 ,501 1 )  ISNRM, (DSN 1 1) , 1  =  1 ,8) 

5011  FORMAT (II ,3X ,8F1 0 .0) 

NUMBER  OF  RPINGS  FOR  THIS  RUN 

READ (5 ,5012)  IRM 

5012  FORMAT (15) 


DESIRED  PARAMETERS  FROM  INPUT  FILE 
1 0  I  SNR- NW~ 0 

READ ( 5 , 5000 )  NF , N , A , B , C , THD , BEAR , DS , DT , R , RMSE , 

#  NS , FOP , PWP , FSP , DSNR , BW, EEES 

5000  FORMAT (I3,2I2,3X,I4,6X,6F10.0,/,3F10.0,4I2,2X,4F10.0,/,2F10.0) 
IF  (EOF  (5)  .NE  .0. )  999, '15 
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15  WR I TE (6 ,6000)  NF , N , A , B ,C , THD , BEAR , DS , DT , R , RMSE , 

M  NS , FOP ,  PWP  ,  FSP ,  DSHR ,  BW,  EEES 

6000  FORMAT (1  HI , 1 X , MS t NQLE  RUN  INFORMATIONM,//3X,MSCHOOL  NUMBER  , 

#  13,212, 3X, MIN, A, B,C  114)  :  AMP  (12)  :  DOPP  II2))M, 

*  /5X.MN  =M,I5,3X,MA  =X , F6 . 1 , 3X , MB  =M , F6 . 1 , 3X ,MC  =M,F6.1, 

M  3X.MASPECT  =M,F5.1 .3X.MBEARING  =M ,  F5 . 1 , /5X , MDS  =M , F7 . 1 , 3X , 

M  MDT  =M,F5.1 ,3X,MR  =M , F7 . 1 , 3X , MMSE  =M , F7 .3 , //3X ,MSAMPLE  M, 

M  MPROCESSI NO  NUMBER  =  M,4I2,3X,M(F0,PW,FS  (12)  :  FLTR  (t2)  :M, 

M  M  SNR  (12)  :  NSEED  l!2))M, 

M  /5X.MF0  =M,F6.1 ,3X,MPW  =M,F7.3,3X, 

M  MFS  =M,F7.3,3X,MD  SNR  =M , F6 . 1 , 3X ,MBW  ~M , F7 . 1 , 3X ,MEEES  sm,F7.3) 
F0=F0PM1000. 

PW=PWP/1000. 

FS=FSPM1000. 

I F (BW.QT , FS/2 . )  GOTO  999 
)  PRINT  CODE 

READ ( 5 ,5005 )  NP 
5005  FORMAT  1311) 


MPREVIOUS  RUN  PARAMETERS  FROM  TAPE  STORAGEM 


(READ  UNTIL  SCHOOL  NUMBERS  MATCH  OR  EOF) 

16  READ ( 1 )  NF1 , THD 1 ,BEAR1 , DS 1 , DTI , R 1 ,RMSE1 ,NS1 ,ASNR1 ,DSNR1 ,BW1 , 

M  L.EEES1 ,VMN,RLTE,RLT1 ,VMX 
IF (EOF (1 ) .NE.O. )  GOTO  999 

WR I TE (6 , 6001 )  NF1 , THD1 .BEAR1 ,0S1 , DTI ,R1 ,RMSE1 ,NS1 , DSNR1 , BW1 , L, 

M  EEES1 

6001  FORMAT! //I X,M INPUT  TAPE  FORMATM,//3X,MSCHOOL  NUMBER  =M, 13,212,3 
M  X , MASPECT  =M,F5. 1 ,3X,MBEARING  =M , F5 . 1 , /5X ,MDS  =M , F7 . 1 , 3X , MDT  =M 
M  , F5 . 1 , 3X ,MR  =M,F7,1 .3X.MMSE  =M , F7 . 3 , //3X .MSAMPLE  PROCESSING  M 
M  MNUMBER  =M , 4 1 2 , 3X ,MD  SNR  =M , F6 . 1 , 3X , MBW  =M , F7 . 1  , /5X ,ML  =M,I3,3X 
M  , MEEES  =M,F7.3) 

I F  I NF1  ID  .EQ.NF(D)  GOTO  18 
READ ( 1 ) 

i F (EOF  II) .NE.O.)  GOTO  999 
READ ( 1 ) 

I F (EOF ( 1 ) .NE.O.)  GOTO  999 
READ  ( 1 ) 

I F (EOF  1 1 ) .NE.O.)  GOTO  999 
READ ( 1 ) 

1 F (EOF  1 1 ) .NE.O.)  GOTO  999 
GOTO  16 


MMMFISH  SCHOOL  SECT  I ONMMM 

18  D=DS-DT  %  D1=DS1-DT1 

SH= . 057  $  SVs.086  $  SHP=5.7  9  SVP=8.6 
GH=SHMTWOPIMFO/CW  9  GV=SVMTWOP IMFO/CW 
WR I TE (6 , 6002) 

6002  F0RMATI//1X.MRUN  PROGRESS IONM , //2X ,MF I SH  COODINATES  AND  M 
M  MCHARACTERISTICS  SECTIONM) 

C  IF  THE  ORIENTATION  IS  THE  SAME,  GOTO  70 

IF (THD1 . EQ.THD.A.BEAR1 . EQ.BEAR.A.D1 . EQ.D.A.R1 .EQ.R)  GOTO  70 
Nw^  -j 

THL= (THD-THD1 ) /RADEQ  9  RMSEL=RMSE 
THSsSINITHD/RADEG)  9  THC=C03 ITHD/RADEG) 

C  SCHOOL  DIMENSIONS  FROM  ELLIPSOID 
THVE:RADEQMATAN2(D,R) 

X1=RMTHC  9  Y1=-RMTHS 
A2=AMM2  9  B2-BMM2 
DHEsRMANGLEIXI ,Y1) 

XI =-XlMl 00 ,  9  Y1=Y1M100. 

DRE“RM1 00 .MANGLE  I Y 1 ,X1) 

DVE=2.MC 
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BR=BEAR/RADEG 

SHS=SH/2.#SIN(BR)  $  SHC=SH/2 .#COS (BR) 

C  IF  NO  CHANGE  IN  THE  ASPECT  ANGLE,  GOTO  20 
IF (THL.EQ.O. I  GOTO  20 
THLS=SINITHL)  S  THLC=COS(THL) 

WRITE (6 ,6005) 

6005  FORMAT (/3X,#NEW  ORIENTATION:  NEW  COORD  I  NATES , DELAYS ,  BEAM  # 

#  ^CORRECTIONS,  AND  MULTIPLE  SCATTERING  LOSSES#) 

READ  ( 1 1 

I F (EOF ( 1 ) . NE . 0 . )  GOTO  999 
READ ( 1 )  X , Y , Z , T , BCL , BCU , RMS , AMP , EFL , EFU , DP 
I F (EOF  1 1 ) . NE . 0 . 1  GOTO  999 
DO  19  1=1 , N 
XL=XII) 

X ( I ) =XL#THLC-Y ( I )#THLS 

19  Yll)=XL#THLS+YII)#THLC 
GOTO  25 

20  WRITE (6 ,6010) 

6010  FORMAT l/3X,#NEW  ORIENTATION:  NEW  DELAYS,  BEAM  CORRECTIONS,  AND# 

#  #  MULTIPLE  SCATTERING  LOSSES#) 

READ ( 1 1 

IFIEOF11) . NE . 0 . )  GOTO  999 
READ (1)  X ,  Y , Z , T , BCL , BCU , RMS , AMP , EFL , EFU , DP 
1 F ( EOF ( 1 ) .NE.O. )  GOTO  999 
25  DO  60  1=1 , N 
XR=R-X ( 1 1 
ZD= (Z ( I ) +D)##2 
XYZ=SQRT (XR##2+Y ( I )##2+ZD) 

#PROPAGATION  DELAYS# 

Til , 1 ) = (SORT (IXR-SHS)##2*IYII)+SHC)##2+ZD) +XYZ) /CW 
Til  ,2  )  =  (SORT  I  IXR+SHS)##2+m  I  )-SHC)##2+ZD)+XYZ)/CW 
Til , 3 ) = ISQRT (XR##2+Y ( I )##2+ (Z ( I ) +D-SV/2 . ) ##2) +XYZ) /CW 
T 1 1 , 4)  =  (SORT  I XR##2+Y ( I ) ##2  + (Z ( 1 1 +D+SV/2 . ) ##2 ) +XYZ) /CW 

^CALCULATION  OF  BEAM  CORRECTIONS# 

IF  ANY  FISH  IS  IN  A  SIDE  LOBE,  A  NEW  VARIABLE  IS  NEEDED 
TO  CHANGE  THE  SIGN  OF  AMP  1 1 )  FOR  L-R.U-D. 

THH= (BR-ATAN2 ( Y 1 1 1 ,XR))##2 
THV= (ATAN2 IZ ( I ) +D , SORT (XR##2  + Y 1 1 )##2) ) )##2 
THVEQS=SQRT(THH#.4386tTHV)#1 0.8071 3 
THHEQS=SQRT ITHH+THV#2 . 28) #7 . 1 62545 
THVEQ=SQRT ( THH# . 1 1 0+THV) #1 0 .80713 
THHEQ=SQRT (THH+THV# . 549 ) #7 . 1 62545 
IFITHVEQS.LT. .00001 .0. THVEQ.LT. .00001)  GOTO  30 
El = . 01 71 89#THVEQS+ . 2 
E2= . 01 71 89#THVEQ+ . 2 

BCL ( I ) = ( (ABSIS IN (THVEQS)/THVEQS) )##E1 + 

BCL (I)  =  1 1 ABS IS  I N ITHVEQS) /THVEQS) )##E1 ♦ 

#  (ABS  CSI N (THVEGi/THVEQ) )##E2-2. )#45 . 

GOTO  35 

30  BCL ( I ) =0 . 

35  IFITHHEQS.lt. .00001 .0. THHEQ.LT. .00001 )  GOTO  40 

El = .0171 89#THHEQS* . 2 
E2=.017189#TKHEQ+.2 

BCU ( 1 )  =  ' (ABS (5 1 N (THHEQS) /THHEQS) )##E1  + 

#  (ABSIS IN (THHEQ) /THHEQ) ) ##E2-2. )#45. 

GOTO  45 

40  BCU(l)=0. 

45  CONTINUE 

IF  FIRST  ORDER  MULTIPLE  SCATTERING  IS  USED, 

THE  NEXT  LINE  IS  CHANGED. 
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RMS ( I ) =RMSEL 
60  CONTINUE 

CALL  STAT(X,N,XM,SX) 

CALL  STATIY.N, YM.SY) 

CALL  STAT ( Z , N , ZM, SZ) 

CALL  STAT (BCL , N , BCLM, SBCL ) 

CALL  STAT (BCU , N , BCUM , SBCU) 

CALL  STAT ( RMS, N.RMSM.SRMS) 

C  SCHOOL  DIMENSIONS  FROM  COORDINATES 
XR-R-XM 

THHC=BEAR-RADEGKATAN21YM,XR) 

ZD=D*ZM 

XYZ=SQRTIXR#)*2  +  YM)*)*2) 

THVC=RADEG#ATAN2 (ZD , XYZI 

DHC=SR5#SY 

DVC=SR5*SZ 

DRC=SR5#SX 

GOTO  80 


NO  CHANGES  IN  ORIETATION 
70  RMSEL=RMSE/RMSE1 
WR I TE 16 ,6015) 

6015  FORMAT  I /3X.HSAME  ORIENTATION)*) 

READ  1 1 )  THVE , DHE , DVE , DRE , THVC , THHC , DHC , DVC , DRC , XM , SX , YM , SY , ZM , SZ , 
«  AMPM , SAMP , BCLM , SBCL , BCUM , SBCU , RMSM , SRMS , EFLM , SEFL , EFUM , SEFU , 
u  n pm  Qnp 

1 F I  EOF  Cl) , NE.O.)  GOTO  999 

READ ( 1 )  X , Y , Z , T , BCL , BCU , RMS , AMP , EFL , EFU , DP 

IF (EOF (1 ) .NE.O.)  GOTO  999 

DO  75  1=1  ,  N 

RMS  1 1 ) =RMS ( I ) XRMSEL 
75  CONTINUE 

IF  NO  CHANGE  IN  AMPLITUDES,  GOTO  100 
80  IFINF1 (2) .EQ.NFI2))  GOTO  100 

^SCATTERING  AMPLITUDE  CALCULATIONS)* 

NW=  1 

WRITE (6 ,6025) 

6025  FORMAT (/3X.XNEW  AMPLITUDES)*) 

C  IF  EQUAL  AMPLITUDES  ARE  DESIRED,  GOTO  100 
IFINFI2) . EQ . 0)  GOTO  90 
C  USE  SEED  ON  FIRST  SCHOOL  ONLY 
IF  C IR.GT. 1 )  GOTO  82 
NR-NF ( 2 ) 

CALL  RANSETINSEED (NR) ) 

82  DO  85  1=1  , N 

FM=RANF (0 . ) ♦ . 5 
RN=RANF (0 . ) 

RN=AMAX1 (RN,. 00001) 

AMP  ( I )  =SQRT  ( -ALOG I  RN )  )*  ( 2 .  #FM)*I*2 ) ) 

85  CONTINUE 
GOTO  96 

C  EQUAL  AMPLITUDES  FOR  ALL  FISH 
90  DO  95  1=1  , N 

AMP ( I ) =1  . 0 

95  CONTINUE 

96  CALL  STAT ( AMP, N, AMPM, SAMP) 

GOTO  105 

100  WRI TE (6 ,6030) 

6030  FORMAT  I /3X,)«SAME  AMPLITUDES)*) 

C  IF  NO  CHANGE  IN  DOPPLER  SHIFTS,  GOTO  125 
105  IF (NF1 (3) .EQ.NFI3) )  GOTO  125 


ooo  oo  o  o  o  ooo 


70 


#D0PPLER  CALCULATIONS# 

WR 1 TE ( 6 , 603  5 ) 

6035  FORMAT (/3X,#NEW  DOPPLER  DATA#) 

IF  NO  DOPPLER  SHIFTS  DESIRED,  GOTO  115 
!F(NF(3) .EQ.O)  GOTO  115 
IFMR.GT.1)  GOTO  107 
NR=NF 13 ) 

CALL  RANSET (NSEED (NR) J 
107  DO  110  1=1 , N 

DP ( I ) = IRANF (0 . )  + . 25 )#FO/375 . 

110  CONTINUE 
GOTO  121 

NO  DOPPLER  SHIFTS  USED 
115  DO  120  1=1, N 
DP  1 1 ) =0 . 

120  CONTINUE 

121  CALL  STAT (DP , N , DPM, SDP) 

GOTO  130 

125  WRI TE (6 , 6040) 

6040  FORMAT (/3X,#SAME  DOPPLER  DATA#) 

IF  NO  CHANGE  IN  ECHO  ENERGIES,  GOTO  140 
130  IF (NW.EQ.O)  GOTO  140 

TOTAL  L-R  AND  U-D  RECIEVED  ENERGIES  (CORRECTED  AMPLITUDES  SQUARED) 
FOR  EACH  FISH 
DO  135  1=1, N 

XYZ=AMP 1 1 ) #RMS ( I ) 

E=BCL( I )/20. 

EFL  ( I )  =  (XYZ#10.##E)##2 
E=BCU( l)/20. 

EFU ( I ) = (XYZ#1 0 ,##E)##2 
135  CONTINUE 

CALL  STAT(EFL,N,EFLM,SEFL) 

CALL  STATIEFU,N,EFUM,SEFU) 

140  CONTINUE 

###ECHO  PROCESSING  SECTION### 

WR  I TE (6 , 6045) 

6045  FORMAT ( //2X , #QUAD  SAMPLE,  FILTER,  AND  NOISE  SECTON#) 

C  IF  THE  SAME  SAMPLES  ARE  USED,  GOTO  200 
IFINS1  (1)  .EQ.NS(D)  GOTO  200 
WRITE (6 ,6050) 

6050  FORMAT l/3X,#NEW  UNFILTERED  QUAD  SAMPLES#) 

C  CALCULATE  THE  NUMBER  OF  SAMPLES  IN  UNFILTERED  ECHO 
VMX=AMAX1 ITI1,1),T(1,2),T(1,3),T(1,4)) 

VMN=AMIN1 (T(1 ,1) ,T(1 ,2) , T 1 1 ,3) , T ( 1 ,4) ) 

DO  145  1=2, N 

VMX1 =AMAX1 (T(I,1),T(I,2),T(I,3),T(I,4)) 

VMN1 = AM  INI (T(I,1),TII,2),T(I,3),T(I,4)) 

VMX=AMAX1 (VMX1 ,VMX) 

VMN=AMIN1 ( VMN1 ,VMN) 

145  CONTINUE 
VMX=VMX+PW 
LMX= INT IFS#VMX) 

LMN= INT(FSnVMN) +1 
L=LMX-LMN+1 
IF(L.GT.IOO)  GOTO  908 
RLT=LMN/FS 
RLT1 =RLT 

C  CALCULATE  EXPECTED  UNFILTERED  SIGNAL  AND  NOISE  ENERGIES 
XI =N#A I  NT (PW#FS+ , 000001 )#.375 
EESL=EFLM#X1 
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EESU=EFUM#X1 

EEN=2.#L 

EES=EEES 

C  IF  NO  EXTERNAL  EXPECTED  ENERGY  (EEES)  IS  SPECIFIED, 

C  THEN  USE  THE  MEAN  OF  THE  L-R  AND  U-D  EXPECTED  ENERGIES. 
IF (EEES.EQ.O. )  EES- (EESL+EESU) /2 . 

CALL  NGAIN(EES,EEN,DSNR,G) 

HLOOP  TO  GENERATE  UNFILTERED  QUAD  SAMPLES# 

GENERATE  SIGNAL  COMPONENTS 
FSL  =  1 . /FS 
DO  147  1=  4 

147  ES ( I ) =  EN 1 1 J  =0 . 

EAL=EAU=DAL=DAU=0. 

DO  167  1=1, L 
DO  150  J  =  1 ,4 

150  QS ( I , J , 1 ) =QS ( I , J ,2) =0 . 

SF 1 1 ) =X1 =Y1 =0 . 

DO  160  J=1 , N 

XR=SQRTf Er L ( J ) ) 

XYZ=SQRT (EFU ( J ) ) 

DO  155  K=1 ,4 

PU=PULSE(RLT-TCJ,K) .J) 

IF (PU.EQ.O. )  GOTO  155 
YR=TWOP  I#  IFO  +  DP  ( J )  I *CT ( J  . K ) 

IFIK.GT.2)  XR=XYZ 
QStl  ,K,1 )=QSI I ,K, I )+XRXPU#COS(YR) 

QS ( I , K , 2) =QS ( I ,K , 2) +XR#PU#S I N I YR) 

155  CONTINUE 

IF (PU.EQ.O. )  GOTO  160 
SFm=SF(l)+1. 

XI =X1 +EFL I J ) 

Y1 =Y1 +EFU I J ) 

160  CONTINUE 

DO  165  J=1 ,4 

ESS ( I  ,  J)=QSM  ,J  ,1)##2+QSM  ,  J  ,2)##2 
165  ES ( J ) =ES I J ) +ESS (I  , J ) 

RLT=RLT+FSL 
EAL=EAL+X1 
EAU=EAU+Y1 
DAL=DAL+X1##2 
DAU=DAU+Y1##2 
SDAII ,1)=X1 

167  SDAII, 2)=Y1 

C  GENERATE  NOISE  COMPONENTS 
IFIIR.GT.1)  GOTO  168 
NR=NS(4) 

CALL  RANSET (NSEED (NR) ) 

168  DO  175  1=1, L 

DO  170  J=1 , 4 
RN=RANF (0 . ) 

RN=AMAX1 (RN, .00001) 

XI =SQRT I -2 ,#ALOG(RN) ) 

Y1 =TWOP l#RANF (0. ) 

QN ( I , J , 1 ) =X1 #COS I Y1 )#G 
QN ( I ,J,2)=X1#SIN(Y1)#G 
XI =QN ( I , J , 1 )##2*QN ( I , J , 2)##2 
ESS  1 1 , J ) =AMAX1 (EGS II , J ) , . 000001 ) 

SSNR 1 1  ,  J )  =1 0 .#AL0G1 0 (ESS  1 1 , )  /XI ) 

170  EN (J ) =EN (J ) +X1 
175  CONTINUE 

RLTE=RLT-FSL 

C  IF  FILTER  IS  DESIRED,  GOTO  180 
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IF (NS (2) . NE . 0)  GOTO  180 
WRITEI6, 6055) 

6055  FORMAT I/3X.XN0  FILTER  APPLIED)*) 

IFINF1 (2) . EQ. -1 )  GOTO  350 
READ ( 1 ) 

I F (EOF ( 1 ) . NE . 0 . )  GOTO  999 
READI1 ) 

IF (EOF (1 ) .NE.O. )  999,350 

XSAMPLE  FILTERING)* 

180  NS1 (2) =0 
READ  ( 1 ) 

I F ( EOF ( 1 ) .NE.O.)  999,201 
C  SAME  CARRIER  FREQUENCY,  SAMPLE  FREQUENCY , AND  PULSE  WIDTH 
C  SO  GET  ENERGIES  FROM  TAPE1  . 

200  READ ( 1 )  ES,EESL,EESU,EN,EEN,EAL, EAU , DAL , DAU 
1 F I  EOF ( 1 ) .NE.O.)  GOTO  999 

C  IF  SAME  FILTERED  ENERGIES  ARE  TO  BE  USED,  GOTO  300. 

IF INS1 (2) . EQ. NS (2) )  GOTO  300 
WRITE (6 ,6060) 

6060  FORMAT (/3X,*APPLY  NEW  FILTER  TO  SAME  SAMPLES*) 

C  TAPE1  SAMPLES  MUST  BE  UNFILTERED  HERE. 

IFINS1 (2) .NE.O)  GOTO  999 
GOTO  202 

201  WR I TE (6 , 6065 ) 

6065  FORMAT (/3X,*APPLY  NEW  FILTER  TO  NEW  SAMPLES*) 

C  CALCULATE  EXPECTED  FILTERED  SIGNAL  AND  NOISE  ENERGIES 
C  CALCULATE  "A”  GIVEN  BW  AND  50  DIVISIONS  IN  HANNING  PULSE. 

202  XR=2 . XCOS ( TWOP I XBWXPW/SO . ) -4 . 

AF=  I -XR-SQRT (XRXX2-4 . ) ) /2 . 

XYZ=1 . -AF 

C  LOOP  TO  CALCULATE  EXPECTED  FILTERED  XMT  PULSE  ENERGY 
EPF=X1=0. 

DO  205  I  si ,50 

XI =AF*X1 +XYZX IS  I N(PI*l/50.)**2) 

EPF=EPF+X1**2 
205  CONTINUE 
1=50 

210  1=1+1 

I F ( I .GT.500)  GOTO  999 
XI =X1 *AF 
WT=X1X*2 
EPF=EPF+WT 

IF (WT.GE . .01XEPF)  GOTO  210 

Y2=EPF/50 . / . 375 

EESL=EESL*Y2 

EESU=EESU*Y2 

EES=EEES 

IF (EEES.EQ.O. )  EES= (EESL+EESU) /2 . 

C  APPLY  FILTER  TO  SIGNAL  AND  NOISE  SAMPLES 
C  CALCULATE  ”A”  GIVEN  BW  AND  F5 
XR=2 . XCOS I TWOP I XBW/FS ) -4 . 

AF= ( -XR-SQRT (XRXX2-4 . ) ) /2 . 

XYZ=1 .-AF 

CUF=SQRT (EEN/2 . /L ) 

AITR=XYZ/  ( 1 .  +AF ) 

'oRsSQRT  (AFR) 

“EN=EEN*AFR 

<. ALL  NGA I N (EES , EEN , DSNR ,G) 

C  >r  TAPE1  SAMPLES  CAN’T  BE  USED,  GOTO  214 
I F (NS1 ( 1 ) . NE . NS ( 1 ) )  GOTO  214 
(vAD ( 1 )  QS,QN,SF,SDA 
iF.rOFM)  .NE.O.)  GOTO  999 
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C  USE  RANF  SEED  ON  FIRST  SCHOOL  ONLY. 

214  IFUR.GT.1)  GOTO  215 
NR-NS ( 4 ) 

CALL  HANSET (NSEED (NR) ) 

CALCULATE  YI-1)  PRE-ECHO  NOISE  FILTER  TERMS 

215  DO  220  J  =  1  ,4 

RN=RANF (0. ) 

RN=AMAX1  (P.N,  .00001) 

XI  = SORT ( -2  .  KALOG ( RN) ) 

Y1=TWOPI#RANFIO.) 

QS(1 ,J,1)=QS(1 ,J,1)#XYZ  $  QSI1 ,J,2)=QS(1 ,J,2)#XYZ 
QN 1 1  , J ,  1 ) sxi  NCOS  1 Y1 ) #GUF#GXAFSR#AF+QN 1 1 , J , 1 ) *XYZ*G 
QN ( 1 , J , 2 ) =X1  #S I N I Y1 ) #GUF#G*AFSR#AF+QN ( 1 , J , 2 ) XXYZXG 
ES(J)=QS(1 , J , 1 )#X2+QS (1 , J , 2) X*2 
ESSII , J)SES(J) 

EN ( J ) =QN ( 1 , J , 1 ) XX2+QN ( 1 , J , 2 )#*2 
ES ( J ) =AMAX1 (ESIJ) , .000001) 

220  SSNR ( 1 , J ) =1 0 . XALOG1 0(ES(J)/EN(J)) 

DO  235  1=2,  L 
DO  230  J  =  1  ,4 
XI =Y1 =0 . 

DO  225  K  =  1  ,2 

QS ( I  ,  J  ,K)  =QS  (1-1  ,  J  ,K)#AF+QS(  I  ,  J  , K ) itXYZ 
ON  1 1 ,J,K)=QN1I-1 ,J,K)*AF+aNI I ,J,K)#XYZXG 
XI =X1 +QS ( I , J ,K)##2 
225  Y1 =Y1 +QNI I , J ,K)#X2 
ESS ( I , J ) =X1 
ES(J)=ES(J)+X1 
EN  1 J ) =EN ( J ) +Y1 
XI =AMAX1 (XI , .000001) 

230  SSNR ( I , J ) =1 0 .HAL0G1 OtXI/YI ) 

235  CONTINUE 

C  LOOP  TO  SAMPLE  FILTER  TAIL 
I  =L 

240  1=1+1 

IFI I .GT.100)  GOTO  999 
SSNRP  =  0 . 

C  GENERATE  MORE  NOISE  SAMPLES 
DO  250  J  =  1 ,4 
RN=RANF (0 . ) 

RN=AMAX1 (RN, .00001) 

XI =SQRT (-2 . HALOG(RN) ) 

XI =SQRT ( -2 .XALOG (RN) ) 

Y1 =TWOP IXRANF (0 . ) 

QN ( I ,J,1)sXlNC0S(Y1) 

QN(  I  ,  J  ,2)=X1)lSIN(Y1 ) 

XI =Y1 =0 . 

DO  245  K  =  1  ,2 

QSU,J,K)=QS(M,J,K)#AF 

QN(  I  ,  J  ,K)=QN(  1-1  ,J  ,K)MAF+QN(  I  ,  J  ,K)xXYZNGJ(GUF 
XI =X1 +QN ( I ,J,K)**2 
245  Y1=Y1+QS( I ,J,K)##2 

XI =AMAX1 (XI , .0000001) 

Y1=AMAX1 ( Y 1 , .000001) 

SSNR ( I , J)=10.#ALOG10(Y1/X1 ) 

250  SSNRP=SSNRP+SSNR( I , J) 

C  IF  THE  AVERAGE  SAMPLE  SNR  IS  LESS  THAN  DSNR/10,  GOTO  260. 
I F (SSNRP . LT . DSNR/2 . 5 )  GOTO  260 
DO  255  J=1 ,4 
ESSII , J)=Y1 
EN ( J ) =EN ( J ) +X1 
255  ESIJ) =ES ( J ) +Y1 

SF ( I ) =SDA (1.1) =SDA ( I , 2) =0 . 

RLTE=RLTE+1 . /FS 
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GOTO  240 
260  L  = I  - 1 
GOTO  350 

C  CHANGE  S/N  RATIO  ONLY 

300  WR I TE (6 ,6070) 

6070  FORMAT l/3X,*SAME  SAMPLES  WITH  NEW  S/N  RATIO*) 

READ ( 1 )  QS ,QN , SF , SDA , SSNR 
IF (EOF (1 ) . NE . 0 . )  GOTO  999 
EES-EEES 

IFIEEES.EQ.O. )  EES= (EESL+EESU) 12  . 

C  IF  DSNR  ON  TAPE1  IS  OK,  GOTO  350, 

IF IDSNR1 .EQ.DSNR)  GOTO  350 

301  CALL  NGA I N (EES , EEN , DSNR , G) 

IF (NS (2) . NE .0)  GOTO  999 
X1=-1O.*ALOG10IG**2) 

DO  310  1=1, L 
DO  305  J  =  1 ,4 

QN(I,J,1)=QN(I,J,1)*G 
QN(I,J,2)=QNII,J,2)*G 
305  SSNR ( I , J ) =  SSNR ( I , J ) +X1 
310  CONTINUE 

DO  315  J  =  1 ,4 
315  ENIJ)=EN(J)*G**2 

C  COMBINED  SIGNAL  PLUS  NOISE  SAMPLES  AND  SCHOOL  CHARACTERISTICS 

350  QLR=QLI=QUR=QUI=0, 

DO  351  1=1,4 

EQ(l)=DQ(l)=0. 

351  CONTINUE 

IFINS13) .EQ.O)  GOTO  355 
WRITE  16,6075) 

6075  FORMAT  1/3 X,*RE SUL TS  CALCULATED  Fi.UM  SIC.'IAL  PLUS  NOISE*) 

DO  354  I =1 , L 
DO  353  J  =  1,4 
X1=0. 

DO  352  K =  1  2 

Ql I , jTk) =QS ( I , J ,K) +QN ( I , J ,K) 

352  XI  =X1 +Q( I , J ,K) **2 
SDQII ,J)=X1 
EQ(J)=EQIJ)+X1 
DQ(J)=DQ(J)+X1**2 

353  CONTINUE 

354  CONTINUE 
GOTO  358 

355  WR I TE (6,6080) 

6080  FORMAT 1/3X,*RESULTS  CALCULATED  FROM  SIGNAL  WITHOUT  NOISE*) 

DO  357  J  =  1  ,4 
EQIJ)=ES(J) 

DO  356  1=1  ,L 

Q( I , J , 1 ) =QS ( I , J , 1 )  $  Q ( I , J , 2 ) =QS (1,3,2) 

SDQII ,J)=ESS(I ,J) 

356  DQ(J)=DQ(J)+ESS( I ,J)**2 

357  CONTINUE 

*L-R  AND  U-D  SAMPLE  CROSSCORRELATIONS* 

358  DO  365  1=1, L 

QLR=QLR+Q( I , 1 , 1 ) *Q( I , 2 , 1 ) +Q( I , 1 , 2)*Q( I ,2 , 2) 

QLI=QLI+Q( I ,1 ,2)*Q( I ,2,1 )-Q( I ,1 ,1)*Q( I ,2,2) 

QUR=QUR+Q( I ,3,1 )*Q( I ,4,1)+Q( I ,3,2)*QI I ,4,2) 

365  QUI=QUI+Q(I , 3 , 2)*Q( I ,4,1) — Q ( I , 3 , 1 )*Q( I ,4,2) 

C  NORMALIZATION 

DALN=EAL**2/DAL 
DAUN=EAU**2/DAU 
CA =  SORT  (QLR**2+QL  1**2 ) 
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CAU=SQRT (QUR##2+QU 1**2) 

CALN=CAL/SQRT(EQ( 1 )#EQ(2) ) 

CAUM=CAU/SQRT(EQ(3)#EQ(4) ) 

DO  390  J=1 ,4 

ES(J)=AMAX1  (ESN)  ,  .000001  ) 

ASNR(J)=10.#ALOG10(ES(J)/EN(J) ) 

DQ( J ) =EQ( J )X#2/DQ( J ) 

ENN(J)=EN(J)/EEN 
ESN  I J) =ES I J) /EES 
EQN(J)=EQ(J)/(EES+EEN) 

IF1NSI3) .EQ.0)  EQN ( J ) =  E S N ( J ) 

C  NORMALIZED  QUAD  SAMPLES 
330  X2 ( J ) =  SQRT ( EQ ( J ) ) 

DO  3S5  1=1  ,L 

SDU()  ,J)=SDQ(  I  ,  J)/EQ(  J) 

Q ( I , J , 1 ) =Q( I , J , 1 ) /X2  ( J ) 

385  Q ( I , J , 2 ) =Q ( I ,J,2)/X2(J) 

390  CONTINUE 
EEQ=EES+EEN 

IF ( NS ( 3 ) .EQ.0)  EEQ=EES 
DO  391  1=1 ,L 

SDANt 1,1) =SDA (1,1) /EAL 

391  SDAN (1,2) =SDA ( I , 2) /EAU 

XSCHOOL  LOCATIONS  AND  DIMENSIONS  FROM  0UAD  SAMPLES# 

THHQ=-ATAN2(QLI ,QLR)#RADEG/GH 
THVQ=-ATAN2 (QU I ,QUR)#RADEG/GV 
IFICALN.LE.1 .)  GOTO  395 
CALP  =  1  . 

GOTO  400 
395  CALP=CALN 

400  IF (CAUN.LE.1 . )  GOTO  405 
CAUP  =  1  . 

GOTO  410 
405  CAUP=0AUN 

410  X1=SR2/GH#SQRT(1 .-CALP) 

Y1=SR2/GV#SQRT(1 .-CAUP) 

DHQM  J=SR5#R#X1 
DVQ ( 1 ) =SR5#R#Y1 
DHQ(2)=X1#RADEG 
DVQ ( 2 ) =Y1#RADEG 
XI =1 800 . /FS 

DRQ ( 1 ) = ( DQ ( 1 ) +DQ (2 ) +DQ ( 3 ) +DQ(4) ) /4 .#X1 
DRQ12) = (DALN+DAUNJ/2 .XXI 

XXXDATA  STORAGE  ON  TAPE  SECTION*## 

WR I TE ( 2 )  NF , THD , BEAR , DS , DT , R , RMSE , NS , ASNR , DSNR , BW, L , EEES , 

*  VMN.RLTE.RLT1 ,VMX 

WRITE (2)  THVE , DHE , DVE , DRE , THVC ,THHC , DHC , DVC , DRC ,XM, SX , YM, SY , ZM, SZ, 

*  AMPM , SAMP , BCLM, SBCL , BCUM , SBCU , RMSM.SRMS , EFLM.SEFL , EFUM, SEFU , 

*  DPM.SDP 

WRITE (2)  X , Y , Z, T , BCL , BCU , RMS, AMP , EFL , EFU , DP 
WRITE (2)  ES,EESL,EESU,EN,EEN,EAL, EAU , DAL , DAU 
WRITE (2)  QS ,QN, SF , SDAN , SSNR , SDQ 

WR I TE ( 4 )  NF , MS , RMSE , THHC , THVC , DHC , DVC , DRC , THHQ , THVQ , DHQ ( 1 ) , DVQ ( 1 ) , 
H  DRQ, EQN ,CAL ,CAU ,CALN,CAUN 
WR ITE (6 , 6085) 

6085  FORMAT(//1X,#DATA  WRITTEN  ON  DISK#) 

XXXPRINT  SECTION*** 

XFISH  LOCATIONS,  DELAY  TIMES  AND  CHARACTERISTICS* 
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IF(NP(1 ) .EQ.O.O. ISNR.GT.O)  GOTO  510 
WR ITE  (6 ,61  00)  (NF(I) ,1=1,3) 

6100  FORMAT  1 1  HI ,  IX ,#F I SH  LOCATIONS,  DELAY  TIMES,  AND  CHARACTERISTICS 

#  #,//3X,#SCHOOL  NUMBER  =#,13,1  OX ,#AMPL I TUDE  SEED  =#,12,1  OX ,#DOPP 

#  LER  SEED  =#, I2,//3X,#N0.#,3X,#X  (M)#,4X,#Y#,7X,#Z#,6X,#TL  (SEC) 

#  #,4X,#TR#,8X,#TU#,8X,#TD#,7X,#AMPL#,8X,#BEAM#,7X,#MULT#,8X,#TOT 

#  AL# , 7X , #DOPP# , /86X , #CORR# , 7X , #SCAT# , 3X , #ENERGY# ,  /) 

DO  500  1=1 ,N 

WRITE (6, 61 05)  I ,X ( I ) , Y II) , Z ( I ) , (T ( I , J ) , J  =  1 ,4) , AMP ( I)  , BCL (I)  , 

#  BCU ( I ) , RMS (l),EFL(l),EFU(l),DP(l) 

6105  FORMAT  MX,  1 4 , 3F8 . 2 , 1 X ,  4F1 0 . 7  ,  F8 . 2 , 2  (2X ,  2F6 . 1  ,  F9 . 2 ) ) 

500  CONTINUE 

WRITE (6, 61 10)  XM, YM , ZM, AMPM , BCLM, BCUM, RMSM, EFLM , EFUM, DPM, SX , SY , SZ , 

#  SAMP , SBCL , SBCU , SRMS , SEFL , SEFU , SDP 

6110  FORMAT l /I X , #MEAN# ,41X,F8.2,2(2X,2F6.1 , F 9 . 2 ) , /2X,#SD  #,3F8.2, 

#  3F8.2,41X,F8.2,2I2X,2F6.1 ,F9.2) ) 

#NORMAL I  ZED  SIGNAL  (OR  SIGNAL  PLUS  NOISE)  SAMPLES  AND  CHARACTERISTICS# 

510  IFINP12) .EQ.O)  GOTO  545 

WR I TE ( 6 , 61 1 5)  VMN.RLTE , RLT1 ,VMX 

6115  FORMAT (1  HI , 1 X ,#NORMAL I  ZED  QUADRATURE  SAMPLES  AND  SAMPLE  # 

#  #CHARACTER I  ST  I CS# , //3X , #T I  ME  OF  FIRST  RETURN  =#,F7.4,10X, 

#  #T I  ME  OF  LAST  SAMPLE  =# , F7 , 4 , /3X ,#T IME  OF  FIRST  SAMPLE  =#, 

#  F7 . 4 , 1 0X ,#END  OF  LAST  RETURN  =#,F7.4) 

IF ( NS ( 2 ) .EQ.O)  GOTO  515 

WR I TE ( 6 , 61 20)  BW 

6120  FORMAT (/3X,#F I LTER  BANDWIDTH  =# , F6 . 1 , 2X , #HZ#) 

515  IF ( NS ( 3 ) .EQ.O)  GOTO  530 

SIGNAL  PLUS  NOISE  QUAD  SAMPLES 
WR I TE (6 , 61 25)  NS ( 4) 

6125  FORMAT ( /3X ,#SAMPLES  OF  SIGNAL  PLUS  NOISE#, 10X,#N0ISE  SEED  =#,12 

#  ,//3X,#NO.#,7X,#LEFT#,12X,#RIGHT#,11X,#UP#,14X,#DOWN#,10X,#SNR# 

#  , 4X ,#F I SH/# , 1 0X ,#DENS I T I ES# , /82X ,#SAMPLE# , 3X ,#F I SH  AMPL#,3X, 

#  #QUAD  SAMPLES#,/) 

DO  525  1=1 ,L 

SSNRP=SDQP=0. 

DO  520  J=1 ,4 

SSNRP=SSNRP+SSNR ( I , J) 

520  SDQP  =  SDQP*SDQM  ,  J) 

SSNRP=SSNRP/4 . 

SDQP  =  SDQP/4 . 

SDAP=(SDAN( I ,1 )+SDAN( I ,2))/2. 

WRITE (6, 61 30)  I , ( (Q( I , J ,K) ,K=1 ,2) , J=1 ,4) , SSNRP , SF ( I ) ,SDAP,SDQP 
6130  FORMAT (IX, 1 4 , 4 (2X , 2F7 . 4) , 4X , F6 . 1 , 3X , F4 . 0 , 2F1 2 .3) 

525  CONTINUE 
GOTO  545 
SNR= INFINITY 
530  WR I TE (6 , 61 35) 

6135  FORMAT (/3X,#SAMPLES  OF  SIGNAL  ONLY# , //3X ,#NO .# , 7X , #LEFT# , 1 2X , 

#  #R IGHT# , 1 1X ,#UP# , 14X ,#DOWN# , 7X ,#F I SH/#, 1 0X ,#DENS I T I ES# , /72X , 

#  #SAMPLE#,3X,#F ISH  AMPL#,3X,#QUAD  SAMPLES#,/) 

DO  540  1=1 ,L 

SDQP=0 . 

DO  535  J  =  1  ,4 

535  SDQP=SDQP+SDQ( I , J)/4. 

SDAP= (SDAN ( I , 1 ) +SDAN ( I , 2) ) /2 . 

WRITE (6, 61 40)  I , ( (Q( I , J , K) ,K=1 ,2) , J=1 ,4) , SF ( I ) , SDAP , SDQP 
6140  FORMAT (IX, I4,4(2X,2F7.4) ,3X,F4.0,2F12.3) 

540  CONTINUE 

#UNN0RMAL I  ZED  SIGNAL  AND  NOISE  COMPONENTS  OF  QUAD  SAMPLES# 


545  IF(NP(3) .EQ.O)  GOTO  555 
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WRITE t6 ,61 45)  NS ( 4) , DSNR 

6145  FORMAT (1  HI .  IX ,#UNN0RMAL I  ZED  SIGNAL  AND  NOISE  COMPONENTS  OF  # 

#  ^QUADRATURE  SAMPLES# , //3X,#NOI SE  SEED  SN.I2.5X, 

#  NDESIRED  S/N  RATIO  =#,F6. 2) 

IFINSI3) .NE.O)  GOTO  546 
WRITEI6.6146) 

6146  FORMAT (/3X,#NO I SE  NOT  INCLUDED  IN  SAMPLE  PROCESSING#) 

546  WRI TE (6 , 61 47) 

6147  FORMAT I/3X ,#NO .# , 1 3X ,#LEFT# , 27X ,#R IGHT# , 26X , 

#  #UP# , 29X ,#DOWN# , /1 1 X ,#S IGNAL# , 9X ,#NOI SE# , 1 1X ,#S IGNAL# , 9X ,#NOI SE 

#  #, 1 IX, #S IGNAL#, 9X, #NOI SE#,1 IX, #S IGNAL#, 9X.NNOI SEN, /) 

DO  550  1  =  1  ,L 

WRITE (6, 61 50)  I , ( ( (QS( I ,J ,K) ,K=1 ,2) , (QNt 1 ,J ,K) ,K=1 ,2) ) ,J=1 ,4) 
6150  FORMAT (IX, 1 3 , 4 (2X , F6 . 3 , F7 .3 , 1 X , F7 . 4 , F8 . 4) ) 

550  CONTINUE 

# INPUT  DATA  AND  COMPUTED  RESULTS# 

555  WRITE (6, 61 55)  NF ( 1 ) , N , SHP , NF 12) ,SVP,NF(3) , FOP, R,PWP,DS,DT, BEAR, 

£  ABC  THD  FSP 

6155  rORMATMHI , IX,# INPUT  DATA  AND  COMPUTED  RESULTS# , //3X , 

#  y’hyUT  DATA# , //4X , NSCHOOL  CHARACTERISTICS#, 30X,#TRANSDUCER  # 

#  >L;,.*,„.\CTERISTICS#,//5X,#SCHOOL  NO.  =#,12X,  13, /5X, 

#  #*\  ■}”  FISH  =# ,  1 0X ,  1 4 , 25X , 

#  #HOJ1‘ .ONTAL  SEPARATION  =#,5X,F4.1,#  CM#,/7X,#AMPLITUDE  SEED  = 

#  # ,7X , 12, 25X ,#VERT ICAL  SEPARATION  =#,7X,F4.1,#  CMN./7X, 

#  NDOPPLER  SEED  =#,9X, I2,25X,#FREQUENCY  =# , 1 5X , F8 . 3  ;  #  KHZ#, 

#  /5X.NRANGE  s#,15X,F7.1 ,#  M# , 1 9X ,#PULSE  LENGTH  s#,14X,F6.3, 

#  #  MSEC#, 

#  /5X, NSCHOOL  DEPTH  s# , 9X , F6 . 1 ,#  M#, 19X,#TRANSDUCER  DEPTH  =1* , 8X 

#  , F6 . 1 ,#  M#, 

#  /5X.NELLIP  SO  I D  PARAMETERS#, 32X.NBEARING  TO  SCHOOL  (+sCCW)  s#, 

#  F5 . 1 ,#  DEG#, 

#  /7X , 1 3HA  (NXNAXIS)  s,9X,F5.1,#  M#,/7X,13HB  (#Y#AX IS)  s,9X, 

#  F5 . 1  #  M# 

#  /7X|#C  (Z  Axis)  =#,10X,F5.1 ,#  M#,/7X,#ASPECT  ANGLE  (+=CCW)  =# 

#  F5 . 1 ,#  DEG#, 

#  //4X.NSAMPLER  CHARACTERISTICS#, //5X.NSAMPLE  FREQUENCY  =#,5X, 

#  F8 .3 ,#  KHZ#) 

IFINSI2) . EQ.O)  GOTO  556 
WRITE (6 .6156)  BW 

6156  FORMAT (5X,#F I LTER  BANDWIDTH  s#,5X,F6.1,N  HZ#) 

GOTO  557 

556  WR I TE (6 , 61 57 ) 

6157  FORMAT (5X,#NO  FILTER#) 

557  IF ( NS ( 3 ) .EQ.O)  GOTO  558 
WR I TE (6 , 61 58)  DSNR,NS(4) 

6158  FORMAT15X.NS/N  RATIO  =# , 1 4X , F4 . 1 ,#  DB#,/5X,#NOISE  SEED  =#, 

#  1 1 X , 14) 

GOTO  559 

558  WR I TE (6 , 61 59) 

6159  FORMAT(5X,#NO  NOISE#) 

C  COMPUTED  SCHOOL  CHARACTERISTICS 

559  WRITE (6 ,6160) 

6160  FORMAT ( //I X ,#COMPUTED  SCHOOL  CHARACTERISTICS  AS  SEEN  FROM  # 

#  NTRANSDUCER# , //3X , MEXPECTED  AND  ACTUAL  COMPUTED  VALUES#) 

I F (NS (3 ) .EQ.O)  GOTO  561 

WR I TE (6 ,61 65)  IESNIJ)  , J  =  1  ,4)  ,EES 

6165  FORMAT (/5X.NENERG I ES  NORMALIZED  TO  THEIR  EXPECTED  VALUES#, 

#  //19X,#LEFT#,5X,#RIGHT#,6X,#UP#,8X,#DOWN#,7X,#EXPECTED#, 

#  //7X ,#S IGNAL# , 4F10.3.F13.3) 

IF (EEES.EQ.O. )  GOTO  560 
WRITE (6 , 61 66) 

6166  FORMAT ( 1 H+ , 3X , # (EXPECTED  ENERGY  NOT  CALCULATED  FROM  THIS  PING)#) 
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560  WR I TE ( 6 , 61 67)  (ENN ( J ) , J= i , 4) , EEN , (EQN (J ) , J*1 , 4) , EEQ, 

X  (ASNRI J) , J  =  1 ,4) ,DSNR 

6167  FORMAT (7X.XNO I SE  X , 4F1 0 . 3 , FI 3 . 3 , /7X .XECHOX , 2X , 4F1 0 . 3 , FI 3 . 3  , 

X  //7X,XSNRX,3X,4F10.3,F13.3) 

GOTO  565 

561  WR I TE (6 , 61 70)  (ESN(J) , J = 1 ,4) ,EES 

6170  FORMAT  I /5X ,xS IGHAL  ENERGIES  NORMALIZED  TO  THEIR  EXPECTED  X, 

X  XVALUEX,//8X,XLEFTX,5XIXRIGHTXI6X,XUPXI8X,XDOWNX17X1XEXPECTEDX, 
X  //2X,4F10.3,F13.3) 

IF (EEES.EQ.O. )  GOTO  565 
WR I TE (6,6171 ) 

6171  FORMAT ( 1 H+ , 3X ,X (EXPECTED  ENERGY  NOT  CALCULATED  FROM  THIS  PING)X) 

C  CORRELATION  AMPL I TUDES 

565  WR I TE 1 6 , 61 75)  CAL,CALN,CAU,CAUN 
6175  FORMAT  I /5X , XCORRELAT I  ON  AMPLITUDESX,//26XlXUNNORMX, 

X  13X .xNORMx , //7X .XLEFT/R IGHTX ,F15.3,F18.6, /7X .XUP/DOWNX , 3X,F1 5 .3 

X  ,F18 .6) 

C  SCHOOL  CHARACTERISTICS 

WRITE  1 6, 61 80)  BEAR , THHC , THHQ , THVE , THVC , THVQ , DUE , DHC , IDHQI I ) , I =1 ,2) 
X  ,DVE,DVC, (DVQ(I) , 1=1 ,2) .DRE.DRC, (DRQ(I) , 1=1 ,2) 

6180  FORMAT 1/3X.XSCHOOL  CHARACTERISTICSX,//5X,XANGLES  TO  SCHOOLx, 

X  /28X,XELLIPSOIDX,8X,XFISH  LOCAT IONSX, 1 0X .XQUAD  SAMPLESX, 

X  //7X ,xHOR I ZONTALX , 3 (FI 6 . 2 , X  DEGX) , /7X ,XVERT ICAL  X , 3 ( FI 6 . 2 , 

x  X  DEGx) , /7X , 

X  /5X.XSCHOOL  DIMENSIONSx,//7X,xHORIZONTALX,2(F16.2,x  M  X), 

X  4X, F8 . 2 ,X  M  X , F8 .2 ,X  DEGX, 

X  /7X.XVERTICAL  X.2IF16.2.XM  X) , 4X , F8 .2 ,X  M  X.F8.2, 

X  X  DEGX./7X, 

X  XDOWN  RANGEX, 

X  2IF16.2.X  M  X) ,4X,2(F8.2,X  M  X) , IX ,X (QUAD , AMP)X , //87X , 

X  1 6HD IM=2xSQRT t 5)XSD) 


IF  ONLY  ONE  SNR  IS  DESIRED  OF  LAST  SNR  HAS  BEEN 
COMPUTED  FOR  THIS  PING,  GOTO  990 
I F C ISNRM.EQ.1 .0. ISNR.EQ. ISNRM)  GOTO  990 
NS (2) =0 
NS (3) =NS (3 ) +1 
I SNRs I SNR+1 
DSNR=DSN ( I  SNR) 

WRITEI6.6185) 

6185  FORMAT (1  HI ,3X,XSAME  SAMPLES  WITH  NEW  S/N  RATIOX) 

GOTO  301 
990  I R- 1 R+ 1 

IFMR.GT.  IRM)  GOTO  999 
GOTO  10 
999  STOP 
END 
C 

Cxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

c 

FUNCTION  PULSE  (T,l) 

c  CALCULATES  AMPLITUDE  OF  SIGNAL  RETURN  FROM  FISH  "I”  AT  TIME  ”T” 

COMMON  /PLSE/PW, FO , DP  1 500) 

PW1 =PWXF0/ CFO+DP  C I ) ) 

IF1T.LT. 0. . 0R.T.GT.PW1 )  1 ,2 

1  PULSE=0. 

GOTO  3 

2  PULSE— (SIN(TX3.141 592654/PW1 ) )XX2 

3  CONTINUE 
RETURN 
END 

C 

C 


oo  o  o  o 
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FUNCTION  ANGLE (XI , Y1 ) 

C  COMPUTES  ANGLE  BETWEEN  TANGENTS  TO  ELLIPSE  FROM  TRANSDUCER 
COMMON  /ANGL/A , A2 , B2 , P I 
SLP=SQRTIB2#X1#*2+A2HY1##2-A2#B2) 

PI1=0. 

SGN1 =  SGN2  =  1 . 

IFIX1 .LT.0.0)  SGN1=-1. 

I F (SGNIXXI -A)  10,30,20 
10  P  1 1 =P  I 
SGN2=-1 .0 

20  SLP1  =  CXtJtYI  -SLP)  / (X1##2~A2) 

SLP2=  (X1*tY1  +SLP) / (X1##2-A2) 

GOTO  40 

30  SLP1  =  IX1#Y1-SLP)/m##2-B2) 

SLP2=  CXI^YI  +SLP)  /  C  Y1J4J42-B2) 

40  ANGLE=SGN2#ABS(ATAN(SLP1 ) -ATAN (SLP2) ) +P 1 1 
RETURN 
END 


SUBROUTINE  STAT I SA , I S , SM, SSD) 

CALCULATES  MEAN  AND  STANDARD  DEVIATION  OF  A  SEQUENCE 
DIMENSION  SA(IS) 

SM=SSD=0. 

DO  10  1=1 , IS 
SM=SM+SA( I ) 

10  CONTINUE 
SM=SM/IS 
DO  20  1=1  ,  IS 

SSD=SSD+(SA( l)-SM)##2 
20  CONTINUE 

SSD=SQRT(SSD/IS) 

RETURN 

END 


SUBROUTINE  NGAIN(ES,EN,DSNR,G) 

C  NOISE  GAIN  TO  ACHIEVE  DESIRED  S/N  RATIO 
SNR=1 0 . XALOG1 0 (ES/EN) 
E=(SNR-DSNR)/20. 

G=10.##E 

EN=ENXG##2 

RETURN 

END 


-END  OF  INFORMAT  I ON- 
? 


APPENDIX  E 


FORMULATION  OF  AMPLITUDE  CORRECTIONS  FOR 
TRANSDUCER  BEAM  PATTERNS 


The  echo  amplitude  of  each  fish  is  partially  determined  by  its 
position  in  the  beam  pattern.  This  appendix  shows  how  these  ampli¬ 
tude  corrections  are  derived  from  the  calibration  curves  of  an 
actual  transducer  used  in  a  split-beam  system.  The  transducer 
is  composed  of  many  elements  which  are  grouped  into  four  quadrants 
as  shown  in  Fig.  13. 


VERTICAL 

CUT 


HORIZONTAL 

CUT 


Figure  13.  Diagram  of  actual  split- 
beam  transducer. 

Quadrants  1  and  2  are  summed  and  quadrants  3  and  4  are  summed 
to  obtain  the  up-down  transducer -half  signals,  while  L  and  4  along 
with  2  and  3  provide  tire  left-right  si  glia  is .  The  transmit  pulse 
originates  from  all  quadrants  simultaneously..  -  -  : 


81 


Horizontal  (90°  cut)  and  vertical  (0°  cut)  calibration  data 
for  the  different  configurations  are  shown  in  Fig.  14.  The 
isolated  points  on  the  calibration  curves  are  the  values  of 
Eq.  22,  which  was  chosen  to  fit  the  curves  and  used  in  the  simu¬ 
lation. 


sin  <j> 


(0.017189  4)  +  0.2) 


(22) 


where 

it  0  XMT  and  left-right 

16.66  0°  cut 

n  0  XMT  and  up-down 

25.13  90°  cut 

(Note:  0  is  in  degrees;  <p  is  in  radians.) 

The  two-dimensional  angular  position  of  each  fish  can  be 
transformed  into  a  one-dimensional  angle  which  has  the  same  beam 
strength  on  the  calibration  curve.  To  accomplish  this  transfor¬ 
mation,  the  beam  patterns  are  assumed  to  be  elliptical  with  con¬ 
stant  eccentricities.  The  eccentricities,  computed  as  the  ratio 
of  the  horizontal  (90°  cut)  to  vertical  (0°  cut)  main  lobe  widths, 
are  1.51,  3.02,  and  0.741,  for  the  transmit,  left-right,  and 
up-down  beams.  These  are  the  values  of  e  in  Eq.  7  which  are 
used  to  obtain  the  equivalent  one-dimensional  angles.  Amplitude 
corrections  are  computed  from  the  equivalent  angles  for  combina¬ 
tions  of  the  transmit  and  left-right  beams  or  the  transmit  and 
up-down  beams  as  described  in  Chapter  3. 


Up-Down  Beam 


Up-Down  Beam 


Figure  14.  Calibration  curves  of  aiv  actual  split-beam  transducer. 

The  isolated  points  are  values  of  the  equations  fitted 
to  the  curves. 
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A  different  set  of  amplitude  corrections  is  used  in  some  of 
the  simulations  where  the  bearing,  ip ,  is  zero,  because  of  errors 
in  calculating  the  beam  patterns.  In  these  simulations,  which  are 
denoted  by  an  asterisk  in  Chapter  4,  the  vertical  angle,  THV,  is 
set  equal  to  zero  for  each  fish.  The  amplitude  corrections  there¬ 
fore  are  a  function  of  only  the  horizontal  angle,  THU.  A  hori¬ 
zontal  or  vertical  equivalent  angle  is  computed  using  THH ,  and 
the  combined  transmit  plus  left-right  and  combined  transmit  and 
up-down  amplitude  corrections,  BCL  and  BCU,  are  calculated.  The 
section  of  the  program  code  that  computes  these  corrections  is  shown 
in  Fig.  15.  The  corresponding  beam  patterns  for  BCL  and  BCU  are 
shown  in  Fig.  16.  The  maximum  difference  in  amplitude  corrections 

for  the  two  sets  of  beam  patterns  ranges  from  0.09  dB  for  the 
small  school  to  2.17  dB  for  the  large  school.  Since  the  effect 


C  CALCULATION  OF  BEAM  CORRECTIONS 

THH«BR-aTaN2 (Y (I )/XK)  i  1HV»ATAN2(0#SQRT(XK**2+Y(l)**2i) 
THV6  Q-SGRT ( (THH/ 1.75 )**2+THV**2 ) *10. 60713 
THHtG*SGRT ( THH++2+ ( 1 . 3  +  THV ) **2 )* 7, 1625464 
IFtTHVEQ.LT.. 00001)  G0TU  30 
E«2.*( .017189+THVEG+.2) 

BCL(I)-( (ABS(SlN(TriVEG)/THVEQ) )**E-1. )*90, 

GOTO  35 
30  BCL ( I ) *0 • 

35  IF(THHEG.LT.. 00001)  GOTO  40 
E*2 • * ( •017189*THHEQ+»2) 

BCUtl)»(  (ABS(SlN(THHEGi/THHEQ)  )**t-l.  )*90. 

GOTO  45 
40  BCUU)*0.~- 
45  CONTINUE 


Figure  15.  Listing  of  the  computer  code  used  in  calculating 
the  second  set  of  beam  losses  for  schools  denoted 
by  an  asterisk  in  Chapter  4. 


Transmit  and  Left-Right 


Transmit  and  Up-Down 


Figure  16.  Combined  transmit  and  receive  beam  patterns  for  the 
second  set  of  amplitude  corrections  (denoted  by  an 
asterisk  in  Chapter  4) . 

of  random  amplitudes  is  found  to  be  insignificant  in  Chapter  4, 
these  differences  in  amplitude  corrections,  although  they  are 
based  on  the  spatial  properties  of  the  fish,  should  not.  cause 
major  changes  in  the  spatial  estimates.  The  dimension  estimates 
may  be  slightly  reduced  because  of  greater  beam  losses  for  fish 
farther  from  the  beam  axis. 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  page  (When  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1  REPORT  NUMBER 

2.  GOVT  ACCESSION  NO. 

3.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (end  Subtitle) 

A  SIMULATION  STUDY  ON  THE  SIZE  AND  ANGULAR 
LOCATION  ESTIMATION  OF  FISH  SCHOOLS  USING 

5.  TYPE  OF  REPORT  4  PERIOD  COVERED 

Technical 

SPATIAL  CORRELATION 

6.  PERFORMING  ORG.  REPORT  NUMBER 

APL-UW  8012 

7.  AUTHORS 

8.  CONTRACT  OR  GRANT  NUMB ER(a) 

Roger  G.  Hendershot 

N00024-78-C-6018 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Applied  Physics  Laboratory 

University  of  Washington 

1013  N.E.  40th,  Seattle,  WA  98105 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  &  WORK  UNIT  NUMBERS 

1  1.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Naval  Sea  Systems  Command 

Department  of  the  Navy 

Washington,  D.C.  20362 

12.  report  date 

August  1980 

13.  NUMBER  OF  PAGES 

90 

14.  MONITORING  AGENCY  NAME  &  ADDRESStH  different  from  Controlling  Office ) 

15.  SECURITY  CLASS,  (ot  th/a  raport) 

Unclassified 

15*.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 

16.  DISTRIBUTION  STATEMENT  (of  this  Report) 

Unlimited 

17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  In  Block  20,  If  different  from  Report) 

18.  SUPPLEMENTARY  NOTES 

Master 's.  Thesis 

r  s  .a  ^ 

19.  'KEY  WORDS  (Continue  on  reverse  aide  If  necessary  and  Identify  by  block  number) 

Acoustic  simulations-  >>  Sonar  false  targets^' 

Acoustic  spatial  correlation  V  Splitjlbeam  sonar,’ • 

High  frequency  acoustics'  Underwater  guidance^  and-eont-ror 

Hydroacoustic  fisheries  assessment,'  \  Volume  reverberation,  l-d"  / 

Reverberation  generator'  fREVGENT  ~  < 

20.  ABSTRACT  (Continue  on  reverse  aide  If'neceaaary  and  Identify  by  block  number) 

A  computer  simulation  based  on  a  point  scatterer  model  has  been 
used  to  investigate  the  problem  of  estimating  the  location  and  dimen¬ 
sions  of  fish  schools  using  split-beam  processing  of  hydroacoustic 
data.  Estimates  are  derived  from  the  temporal  crosscorrelation  of 
signals  received  by  transducer  halves,  the  peak  correlation  amplitude 
being  related  to  the  variance  of  the  scatterer  locations  (the  fish  . 

school's  width)  and  the  time  shift  of  the  peak  being  related'  to  the 

fa 


on  F0HM 

WW  I  JAN  73 


1473  EDITION  OF  I  NOV  65  IS  OBSOLETE 
S/N  0102-LF  014-6601 


0 


UNCLASSIFIED 


SECURITY*  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Entered) 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  P AGE(W7ien  Dal t  Entered) _ 

• -k  mean  of  the  scatterer  locations  (the  *ish  school's  angular  position). 

These  estimates  are  formed  with  vertically  and  horizontally  separated 
transducer  halves  and  combined  with  echo  arrival  time  to  provide  three- 
dimensional  position  and  size. 

It  is  demonstrated  that  fluctuations  in  the  estimated  angular  size 
decrease  as  the  ratio  of  pulse  length  to  echo  length  decreases. 

The  effect  of  additive  noise,  uncorrelated  between  acoustic  channels, 
has  also  been  studied.  Noise  causes  .dth  estimates  to  be  biased  toward 
larger  values,  owing  to  decorrelation  between  the  acoustic  channels.  The 


location  estimates,  however,  remain 
down  to  approximately  -5  dB.  (<s 


unbiased  for  signal-to-noise  ratios 


_ _ :  _UfiCLASS  T  FILE  El - - - 

SECURITY  CLASSIFICATION  OF'YnTS  PA3E0*7«rt  Peie-SHilirstU 


